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An introductory lectures on cosmology at ITEP Winter School for students 
specializing in particle physics are presented. Many important subjects are not 
covered because of lack of time and space but hopefully the lectures may serve as a 
' starting point for further studies. 

00 ' 

O ■ 1 Introduction 

Modern cosmology is vast interdisciplinary science and it is impossible to cover it in any 
considerable detail in five hours allocated to me at this School. The task is even more 
difficult because of different background and level of the participants. Planning these 
lectures, I have prepared the following short list of subjects, which is surely will be made 
much shorter at this lecture course, but hopefully it may be useful for the students who 
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would like to continue studying this field. So the idealistic content could be the following: 

1. A little about general relativity and its role in cosmology. 

2. Four basic cosmological equations and expansion regimes. 

3. Universe today and in the past. 

4. Kinetics in hot expanding world and freezing of species. 

5. Inflation: kinematics, models, universe heating, and generation of density pertur- 
bations and gravitational waves. 
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6. Big bang nucleosynthesis. 

7. Field theory at non-zero temperature and cosmological phase transitions. 

8. Baryogenesis and cosmological antimatter. 

9. Neutrino in cosmology (bounds on mass, oscillations, magnetic moment, and anoma- 
lous interactions. 

10. Dark matter and large scale structure (LSS). 

11. Vacuum and dark energies. 

12. Cosmic microwave radiation (CMB) and cosmological parameters. 

In reality about a half of this plan was fulfilled. At least this lectures could be helpful for 
a first aquaintance with cosmology and as starting point for deeper studies. 

We will start from some non-technical introduction to General Relativity and relations 
between the latter and cosmology, sec. EJ Next we will derive the basic cosmological 
equations ina rather naive way studying motion of non-relativistic test body in spherically 
symmetric gravitational field, sec. El There we also talk about realistic regimes of the 
universe expansion and basic cosmological paramters. In the next section, H] the universe 
history is very briefly presented. Section is dedicated to thermodynamics and kinetics 
in the early universe. Section [6] is dedicated to freezing of species and cosmological limit 
on neutrino mass. Big bang nucleosynthesis is presented in sec. [7J In section [8] the role of 
neutrinos in BBN is described. Neutrino oscillations in the early universe are considered 
in sec. HJ In section [TU] inlationary cosmology is discussed, and the last section [TT] is 
dedicated to cosmological baryogenesis. 

2 Gravity and cosmology 

Two simple observations that the sky is dark at night and that there are shining stars 
lead to the conclusion that the universe is finite in space and time. The first one is the 
well known Olbers' paradox, based on the estimate of the sky luminosity, which in infinite 
homogeneous static universe must be infinitely high. Shining stars should exhaust their 
fuel in finite time and thus cannot exist in the infinitely old universe - thermal death of 
the universe. General relativity (GR) successfully hit both targets leading to the notion 
of expanding universe of finite age, but created instead its own very interesting problems 
which we discuss in what follows. 

Newtonian theory of gravity has an evident shortcoming that it has action-at-a- 
distance property. In other words, gravitation acts instantaneously, at any distance. On 
the other hand, in the spirit of contemporary wisdom interactions are always mediated 
by some bosonic fields and are relativistically invariant. If we wished today to generalise 
Newtonian theory of gravity to relativistic theory we could take, a priori as a mediator 
of interactions scalar, vector, or tensor intermediate bosons, confining ourselves to lower 
spins. 
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Since we know that gravity operates at astronomically large distances, the mass of the 
intermediate boson should be zero or very small. Indeed, massless bosons create static 
Coulomb type potential, U ~ 1/r, while massive bosons lead to exponentially cut-off 
Yukawa potential, U ~ exp(— mr)/r. 

Interactions mediated by vector field are odd with respect to charge parity transfor- 
mation, C-transformation, and as one can see from the vector boson propagator, such 
interactions induce matter-antimatter attraction and matter-matter repulsion, recall elec- 
tromagnetic interactions. Hence vector field cannot mediate attractive gravitational force. 

Scalar and tensor mediators lead to attraction of matter-matter and matter-antimatter 
and both are a priori allowed. According to non-relativistic Newtonian theory the source 
of gravity is mass. Possible relativistic generalisation for scalars should be a scalar quantity 
coinciding in non-relativistic limit with mass. The only known such source is the trace 
of the energy- momentum tensor of matter, Tjf. The relativistic equation of motion for 
scalar gravity should have the form: 

«9 2 $ = 8ttG n T^, (1) 

where Gn is the Newtonian gravitational coupling constant. Such theory is rejected by the 
observed light bending in gravitational field, since for photons: = 0. A small admixture 
of scalar gravity to tensor one, i.e. Brans-Dicke theory [I], is allowed. 

There remains massless tensor theory with the source which may be only the energy- 
momentum tensor of matter, T^ v . In first approximation the equation of motion takes 
the form: 

d 2 V = 8ttG n T^. (2) 

This equation is valid in the weak field approximation because the energy-momentum of 
hfj, u itself should be included to ensure conservation of the total enegy-momentum. 

Massless particles, as e.g. gravitons, must interact with a conserved source. Oth- 
erwise theory becomes infrared pathological. The energy-momentum tensor of matter 
is conserved only if the energy transfer to gravitational field is neglected. Taking into 
account energy leak into gravity leads to non-linear equations of motion and allows to 
reconstruct GR order by order. For a discussion of this approach see papers [2] . 

Historically Einstein did not start from field theoretical approach but formulated gen- 
eral relativity in an elegant and economical way as geometrical theory postulating that 
matter makes space-time curved and that the motion of matter in gravitational field is 
simply free fall along geodesies of this curved manifold. This construction is heavily based 
on the universality of gravitational action on all types of matter - the famous equivalence 
principle, probably first formulated by Galileo Galilei. The least action principle for GR 
was formulated by Hilbert with the action given by 

A = — V- / d 4 x^gR + A m , (3) 

where R is the curvature scalar of four dimensional space-time and and A m is the matter 
action, written in arbitrary curved coordinates. Gravitational field is identified with the 
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metric tensor, g^, of the curved space-time. The curvature is created by matter through 
equations of motion: 

Rfiv — ~^9iivR = SttGnT^, (4) 

where R^ u is the Ricci tensor. There is no space here to stop on technicalities of Riemann 
geometry. A good introduction can be found e.g. in book [3] where one can find definition 
and properties of the Christoffel symbols, r° , Riemann tensor, R^ av p, Ricci tensor, R^ u = 
g^Rfuxup, scalar curvature, R = g^R^, covariant derivatives in curved space-time, D^, 
etc. 

The source of gravity is the energy-momentum tensor of matter taken in this curved 
space-time: 

T, u = 2SA m /Sg^. (5) 

The impact of gravity on matter is included into T^ v due to its dependence on metric and 
in some more complicated cases on the curvature tensors. Let us repeat that the motion 
of matter in the gravitational field is simply the free fall, i.e. motion along geodesies. 

Classical tensor theory of gravity agrees with all available data and is a self-consistent, 
very beautiful and economic theory. It is essentially based on one principle of general co- 
variance, which is a generalisation of Galilei principle of relativity to arbitrary coordinate 
frames. Invariance with respect to general coordinate transformation (which is called gen- 
eral covariance) is a natural framework which ensures vanishing of the graviton mass, m g . 
Even if the underlying classical theory is postulated to be massless, quantum corrections 
should generally induce non-zero mass if they are not prevented from that by some sym- 
metry principle. This is another advantage of tensor gravity with respect to scalar one for 
which no principle which forbids non-zero mass is known. Though quantum gravity is not 
yet understood, it is natural to expect that quantum corrections should induce m g ^ in 
absence of general covariance. 

An important property of equations of motion (jl]) is that their right hand side is 
covariantly conserved: 




Accordingly the energy-momentum tensor must be conserved too: 

D„T^ = 0. (7) 

Here, is covariant derivative, as we have already mentioned. To those not familiar 
with Riemann geometry it may be instructive to mention that covariant derivative appears 
when one differentiates in curved coordinate system, e.g. in spherical one, even in flat 
space-time. From another point of view, covariant derivative in curved space-time, which, 
e.g. is acting on vector field, looks as 

DpV v = <^K - r« V a . (8) 
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It is similar to covariant derivative in gauge theories, because the latter includes gauge 
field, analogous to T^ v . 

According to the Noether theorem, the conservation of follows from the least 
action principle if the matter action is invariant with respect to general coordinate trans- 
formation. So the gravitational (Hilbert) part of the action and the matter part lead to 
self-consistent equations of motion only if general covariance is maintained. 

There is a deep analogy between the Einstein gravity and Maxwell electrodynamics. 
The Maxwell equations have the form: 

Q^pr* = A7T j» (9) 

Owing to anti-symmetry of F^, the l.h.s. is automatically conserved: 

d IM d v F'» EE 0, (10) 
so the current must be conserved too: 

V = 0- (H) 

These two conditions are consistent due to gauge invariance of the total electromagnetic 
action with matter included. 

Einstein was the first who decided to apply GR equations (j4j) to cosmology in 1918 
and was very much disappointed to find that the equations do not have static solutions. 
So an advantage of GR was erroneously taken as a shortcoming. Only after the Friedman 
solution in 1922 which predicted the cosmological expansion [I] and the Hubble discovery 
of the latter in 1929 [5], the idea that our world is not stationary and may have a finite 
life-time was established. 

The distribution of matter in the universe is assumed to be homogeneous and isotropic, 
at least in the early stage, as indicated by isotropy of cosmic microwave background 
radiation (CMB), and even now at large scales. Correspondingly the metric can be taken 
as homogeneous and isotropic one (FRW metric [6]): 

ds 2 = dt 2 - a 2 (t) [f(r)dr 2 + r 2 dVt\ , (12) 

where the function f(r) describes 3D space of constant curvature, f(r) = 1/(1 — kr 2 ). 

The evolution of the scale factor a(t), i.e. the expansion law, is determined by the 
Friedman equations, which follow from the general GR ones for the FRW anzats. The 
derivation is straightforward but quite tedious. We will derive them in the next section 
in very simple but not rigourous way. The derivation may be taken as a mnemonic rule 
to recall the equations in one-two minutes. 

3 Cosmological expansion 
3.1 Basic cosmological equations 

Here we will present an oversimplified derivation of the Friedman equations. Though 
the arguments are subject to criticism, the final results are correct. Let us consider a 
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test body on the surface of homogeneous sphere with radius a(t) and the energy den- 
sity p. The energy conservation condition for the non-relativistic test particle reads 
v 2 /2 = G^Mja + const where M = 4ira 3 p/3. It can be rewritten as: 

\a ) 3 a 2 

This is one of the main cosmological equations. Here is the famous Hubble expansion law, 
that the object situated at distance d runs away from us with velocity proportional to 
the distance, v = Hd. Notice that for d > l/H it runs away with superluminous velocity. 
Superluminous velocities of distant objects are allowed by GR but locally velocities must 
be always smaller or equal to the speed of light. 

Another equation follows from the energy balance of the medium inside the sphere: 
dE = —PdV where E = pV and dE = Vdp + 3(da/a)Vp. Hence: 

p + 3H(p + P) = 0. (14) 

This equation is simply the law of covariant energy-momentum conservation ([7]) in metric 

(USD. 

Problem 1. Derive from eqs. (|T3l) and (TBI) the law for the acceleration of the test 
body: 

a 4n Gn 



a 



(p + 3P) (15) 



A striking feature of equation ( 1151) is that not only energy but also pressure gravitates. 
It is always assumed in canonical theory that p is positive, though pressure may be 
negative. Thus if p+3P < 0, the cosmological expansion would proceed with acceleration, 
i.e. antigravity may operate in cosmological scales. In other words, negative pressure is 
the source of the cosmological expansion. Life is possible only because of that. We believe 
that the universe was in such anti-gravitating state at the very beginning, during the so 
called inflationary stage (see below). Surprisingly it was established during the last decade 
that at the present time the expansion is also accelerating. There existed a simple analogy 
between the universe expansion and the motion of a stone thrown up from the Earth with 
some initial velocity Vq. The speed of the stone drops down with time and it either will 
return back to the Earth, if vq is smaller than a certain value, V\. If vq > v\ the stone 
will never come back. In the last case the stone will either come to infinity with non-zero 
speed or with the vanishing one. All three such regimes could exist in cosmology and 
they were believed to be realised, depending upon the initial expansion velocity. The 
first regime corresponds to the closed universe with p > p c , where p c is the critical or 
closure energy density, see below. The expansion in this case will ultimately turn into 
contraction. The other two regimes correspond to open universe which was expected to 
expand forever. The third one with zero velocity at infinity corresponds to spatially flat 
universe, with k = in eq. f)13p . Now the picture is very much different. Imagine that 
you have thrown a stone from the Earth and first the stone moves with normal negative 
acceleration and after a while starts to move faster and faster as if it has a rocket engine. 
This is exactly what we see in the sky now. It means, in particular, that the spatially 
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closed universe may expand forever. This accelerated cosmological expansion, induced by 
antigravity at cosmological scale is prescribed to existence of mysterious dark energy. It 
is one the greatest unsolved problems in modern fundamental physics. 

Problem 2. We seemingly started from the Newtonian theory but came to the con- 
clusion that pressure gravitates which is not the true in Newtonian case. Where is the 
deviation from Newton? 

Problem 3. Prove that for positive definite energy density, p > 0, any object of finite size 
creates an attractive gravitational force, even if inside such an object pressure may be 
arbitrary negative. 

Problem 4- Prove that any finite object with positive energy density gravitates, so anti- 
gravitational action of pressure can manifest itself only in infinitely large objects. 

Except for equations which determine the law of the cosmological expansion, we need 
an equation which governs particle propagation in FRW metric, i.e. geodesic equation. 
The latter can be written as: 

dV a 

-T^V^V" + curvature term, (16) 



ds 



where V a = dx a /ds and the curvature term is absent in spatially flat universe, when 
k = 0. In what follows we will consider only this case, moreover, the effects of curvature 
are typically small. 

To solve this equation one needs first to calculate the Christoffel symbols for metric 
(j!2p . In 3D flat space they have very simple form: 

T) t = H5), Y% = Ha 2 5 ir (17) 

All other are zero. After that the geodesic equation takes a very simple form: 

p = -Hp (18) 

with an evident solution p ~ l/a(£) which describes red-shifting of momentum of a free 
particle moving in FRW background. In derivation of this result one has to pay attention 
that physical momentum is defined with respect to physical length dl = a(t)dx. 

We can simply derived the same result taking into account the Doppler red-shift of 
the momentum of a free particle induced by the cosmological expansion. Let us take two 
points A and B separated by distance dl. The relative velocity of these two points due to 
expansion is U — Hdl. The Doppler shift of the momentum of the particle moving from 
A to B with velocity v = dl/dt is 

dp = -UE = —HEdl. (19) 

Thus 

p = -HEdl/dt = -Hp. (20) 
Let introduce at this stage the notion of the cosmological red-shift: 

z = a(tu)/a(t) - 1, (21) 
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where t u is the universe age and a(%) is the value of the scale factor today. So the mo- 
mentum of a free particles drops down in the course of expansion as p ~ 1/a ~ 1/(2 + 1). 

Equations (I13H15P are the basic cosmological equations for three unknowns, a, p, and 
P. However, there are only two independent equations. So it is necessary to have one 
more equation describing properties of matter i.e. the equation of state (e.o.s.): P = P(p). 
Usually this equation is parametrized in the simple linear form: 

P = wp. (22) 

The parameter w determines matter properties. For non-relativistic matter pressure is 
negligibly small in comparison with p and in a good approximation we can take w — 0. 
For relativistic matter w = 1/3. There is also one more type of matter (or vacuum) known 
to exist in the universe, for which w = —1. 

However, sometimes the equation of state does not exist but the necessary additional 
relation (not e.o.s.) can be derived from the equations of motion. E.g. for a scalar field: 

D 2 <t> + U'((f>) = (23) 

and one can calculate T M „ and find p and P but P ^ P(p)- 

Problem 5. Calculate T^ v {(j)), p, and P for homogeneous field (f)(t). 



3.2 Expansion regimes 

Here we will present solutions of the cosmological equations for several special cases which 
were/are realised in the universe at different stages of her evolution. We always assume 
that the three dimensional space is flat, i.e. k = 0. As we see below it is true during 
practically all life-time of the universe. 

Let us first consider non-relativistic matter with equation of state P = 0. According 
to eq. ( 1141) the evolution of the energy density is given by: 

p = -3Hp (24) 

and thus p ~ 1/a 3 . The result is evident, it is simply dilution of the number density of 
massive particle at rest. 

The time dependence of the cosmological scale factor, is determined by eq. (TTSI : 

a/a ~ y/p (25) 

and thus in non-relativistic regime a ~ t 2 / 3 and H = 2/3*. 

For relativistic matter the equation of state is P = p/3 and correspondingly 

p = -AH p. (26) 

Thus p drops as p ~ 1 /a 4 and the scale factor rises as a (t) ~ t 1/2 , which means that 
H = 1/21 

The energy density of relativistic particles drops one power of a faster than that of 
non-relativistic ones due to dilution of their number density as volume, 1/a 3 , and red-shift 
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of the particle momentum. That's why relativistic matter dominated in the early universe, 
while at a later stage non-relativistic matter took over. Until last years of the XX century 
it was believed that the universe today is dominated by non-relativistic matter but then 
it was established that the dominant matter is the so called dark energy with equation of 
state close to the vacuum one. 

In the vacuum(-like) regime the energy-momentum tensor is proportional to the metric 
tensor which is the only invariant tensor: 

Tfiy Pvac 9^u- (^7) 

Hence P vac = —p vac and vacuum energy density remains constant in the course of the 
cosmological expansion: p = —3H(p + P) =0. The scale factor in this case rises expo- 
nentially a ~ exp(Ht). 

According to our understanding, all visible universe originated from microscopically 
small volume with negligible amount of matter by exponential expansion with practically 
constant p, see below. 

It is interesting to calculate the causality distance as a function of time, which is equal 
to the light path propagating from an initial moment t\ — to final moment £2- This 
distance can be found from the light geodesic equation: dt 2 — a 2 (t)dr 2 = 0. It can be 
easily integrated to give: 

l y = a(t)J^ y (28) 

Thus for relativistic regime, Z 7 = 2t, for non-relativistic one it is Z 7 = 3t, and for expo- 
nential De Sitter (inflationary stage): Z 7 = H _1 [exp(ift) — 1]. 

As we see, cosmological equations do not have stationary solutions, at least for the 
examples taken. In the case of positive space curvature, i.e. k > 0, and normal matter 
with p ~ l/a", n=3,4, the expansion will ultimately change into contraction. If however, 
p > k/a 2 , e.g. if p is vacuum energy, the expansion may last forever for any k. 

Note, that if the cosmological energy density is dominated by the normal matter the 
Hubble parameter drops down as H ~ 1/t, where t is the universe age. If the universe is 
dominated by vacuum(-like) energy, the Hubble parameter remains constant, 

Problem 6. Find p(a) and a(t) for general linear equation of state, P = wp with 
arbitrary w. Study the case of w < — 1. 

3.3 Cosmological parameters 

Before proceeding further let us say a few words about the natural system of units which 
is used throughout all these lectures. We take speed of light, reduced Planck constant, 
and Boltzmann constant all equal to unity, c = h/2ir = k = 1. All dimensional quantities 
have dimension of length, or time, or (inverse) mass or energy - all the same. For ex- 
ample the Newtonian gravitational constant has dimension of inverse mass, Gn = 1/Mpf, 
M Pl = 1.221 -10 19 GeV = 2.176 • 1(T 5 £; m p = 938MeV = 1.67- 10" 24 g; 
1 GeV" 1 = 1.97 • 10" 14 cm = 0.66 • 1(T 24 , s; 1 eV = 1.16 ■ lO 4 ^ . 
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Homogeneous cosmology is described in terms of the Hubble parameter H = a/a, 
which characterises the universe expansion rate, by the critical or closure energy density 

p c = 3H 2 MjV 8tt (29) 

and by the dimensionless parameter Qj = Pj/p c , which measures the relative contribution 
of the energy density of species of type j into the total energy density of the universe. 
Clearly for spatially flat universe the total energy density is equal to the critical one: 
Qtot = 1, if k = 0. It remains constant in the course of the universe expansion. 
If k 0, then from eq. ( TT31) follows that Q evolves with time as 

1 \ p «o 



ft(a) 



(30) 



f2 J pa 2 

where the index sub-0 denotes the present day values of the corresponding quantities. 

For normal matter pa 2 — > if a — > oo and Q runs away from unity: Q(a) — > if Qq < 1 
and fi(a) — > oo if fio > 1- On the other hand, fi(a) — > 1 when pa 2 — > oo, e.g. for vacuum 
energy at expansion for arbitrary initial Q. 

The present day value of the H characterises by dimensionless parameter h as 

H = 100 h km/sec/Mps, (31) 

where h = 0.73 ± 0.05. The inverse quantity H" 1 = 9.8 Gyr/h ~ 13.4 Gyr is approxi- 
mately equal to the universe age. 

An exact expression for the universe age through the present day values of the Hubble 
parameter and relative energy densities of different forms of matter can be obtained by 
integration of the equation 

a= [87rpG N a 2 /3-k] 1/2 (32) 
After simple algebra one finds: 

1 f 1 dx 
t "=H I = = ( 33 ) 



1-^ + ^ + ^+a; 2 ^ 



where Q t is the total Q and £l rm ,v are respectively contributions from relativistic and 
non-relativistic matter and from vacuum energy. All the quantities in this equations are 
the present day ones; we skipped the sub-index here. 

Problem 7. Derive eq. ( 1331) . Find ty for Q t = = 0, 0.3, 1. Find tjj for Q t = 1, 
Q m = 0.3, and Q v = 0.7. 

The "measured" value of the universe age lies in the interval 

tu = 12 - 15 Gyr, (34) 

found from the ages of old stellar clusters and nuclear chronology. 
The present day value of the critical energy density is : 

= ZH 2 m 2 Pl = Lgg 1() _ 29/i2 g = 1Q 5 ^keV ^ 1() _ 47/i2 Gey4 
8ir era 6 cm 6 

It corresponds approximately to 10 protons per m 3 , but the dominant matter is not the 
baryonic one and in reality there are about 0.5 protons per cubic meter. 
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3.4 Matter inventory 

The relative contributions of different forms of matter into the total energy density are 
obtained from different independent astronomical observations. Here we only present 
their numerical values. For discussion of their measurements in more detail a few extra 
lectures are necessary. 

The total cosmological energy density is very close to the critical one: Q tot = 1 ± 0.02 
as found from the position of the first peak of the angular spectrum of CMBR and the 
large scale structure (LSS) of the universe. 

The usual baryonic matter makes quite small contribution: Qb = 0.044 ± 0.004 as 
found from the heights of the peaks in angular fluctuations of CMB, from produciton of 
light elements at BBN, and from the onset of structure formation with small 5T/T. 

Approximately five time more than baryons is brought by the so called dark matter. 
It is invisible matter with presumably normal gravitational interactions. As is found 
from galactic rotation curves, gravitational lensing, equilibrium of hot gas in rich galactic 
clusters, cluster evolution, and LSS: Qdm ~ 0.22 ± 0.04. 

The rest, Qde ~ 0.76,, is carried by some mysterious substance, which is uniformly 
distributed in the universe induces accelerated cosmological expansion. Its equation of 
state is close to the vacuum one, i.e. w ~ — 1. The existence and the properties of dark 
energy was deduced from dimming of high-z supernovae, LSS, CMB spectrum, and the 
universe age. 

I would like to stress that the different pieces of data and their interpretation are in- 
dependent. It minimises the probability of a possible interpretation error. The numerical 
values obtained in different type measurements are pretty close to each other. 

4 Brief cosmological history 

Universe history can be separated into several epochs, some of them are described by 
established well known physics verified by experiment, some are based on hypothetical 
physics beyond the standard model, and some (little?) are absolutely dark. 

1. Beginning, unknown. Quantum gravity, quantum space-time? Maybe time did not 
exist? It is so called pre-inflationary cosmology. 

2. Inflation, i.e. epoch of exponential expansion of the universe. It is practically 
"experimental" fact.. 

3. End of inflation, particle production. At that period dark expanding "emptiness" 
filled by scalar (or some other) field, inflaton, exploded, creating light and other 
elementary particles. 

4. Baryogenesis. At that time an excess of matter and antimatter in the universe (or 
vice versa) was created. 

5. Thermally equilibrium universe, adiabatically cooled down. Presumably during this 
epoch several phase transitions took place leading to breaking of grand unified sym- 
metry (GUT), electroweak (EW) symmetry, supersymmetry (SUSY), phase transi- 
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tion from free quark-gluon phase to confinement phase in quantum chromodynamics 
(QCD), etc. with possible formation of topological defects and non-topological soli- 
tons. At the phase transitions adiabaticity of expansion could be broken. 

6. Decoupling of neutrinos from electromagnetic part of the cosmo logical plasma. It 
took place when the universe was about 1 second old at T ~ 1 MeV. 

7. Big bang nucleosynthesis (BBN), which proceeded in the time interval from 1 s to 
~ 200 s, and T = 1 - 0.07 MeV. At that time light elements, 2 H, 3 He, 4 ife, and 
7 Li were formed. Theory is in a good agreement with observations. A different 
mechanisms for creation of light elements is unknown. It makes BBN one of the 
cornerstones of the standard cosmological model (SCM). 

8. Onset of structure formation which started when the dominating cosmological mat- 
ter turned from relativistic into non-relativistic one. It took place at the red-shift 



9. Hydrogen recombination, at z ~ 10 3 or T ~ 0.2 eV. At that time cosmic microwave 
radiation (CMB) decoupled from matter and after that it propagated practically 
freely in the universe. After decoupling of matter and radiation baryons begun to 
fall into already evolved seeds of structures created by dark matter (DM). 

10. Formation of first stars and reionization of the universe. 

11. Present time, tjj — 12 — 15 Gyr. 

5 Hot equilibrium epoch 

Usually a system comes to the state of thermal equilibrium after sufficiently long time. 
Paradoxically, in cosmology equilibrium is reached in the early universe when time is short 
but temperature is high and the reaction rates T exceed the cosmological expansion rate, 
H = a/a: 



This condition is fulfilled at high temperatures but bounded by T < a n Mpi, where a is the 
generic value of the coupling constant and n = 1,2 for decays and reactions respectively. 
At lower T the equilibrium is broken due to Boltzmann suppression of the participating 
particles. 

In equilibrium particle distribution functions are determined by two parameters only, 
by the temperature, T and, if the particles do not coincide with antiparticles, by their 
chemical potential, /i. The equilibrium distributions have the well known form: 



cq 



10 4 , T ~ eV. 



r = an ~ a 2 T > H ~ T 2 /Mj 



(36) 




1 



(37) 



exp [{E - y)/T\ ± 1 ' 
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where E = \Jp 2 + m 2 . Equilibrium with respect to the reaction a\ + a 2 + a 3 . . . <-> 61 + 6 2 + 
imposes the following condition on the chemical potential of the participating particles: 



(38) 



Chemical potentials are necessary to introduce in charge asymmetric case to describe 
inequality between number densities of particles and antiparticles: n 7^ n. It is assumed 
usually that in cosmological situation chemical potentials are very small, as follows from 
the observed baryon asymmetry of the universe. However, large lepton asymmetry is not 
excluded and so it may be interesting to consider non-negligible chemical potentials. One 
should keep in mind however, that chemical potential of bosons cannot be arbitrarily 
large. It is bounded by the mass of the boson, \l < m, while for fermions there is no 
upper limit on \i. 

What happens if charge asymmetry in bosonic sector, i.e. in — n) is so large that 
fi = m is not sufficient to realise that? In this case the equilibrium distribution function 
acquires an additional term: 



i.e. Bose condensate forms. Notice that equilibrium distributions are always determined 
by two parameters: T and \i if \x < m or T and C if /i is fixed by the maximally allowed 
value, n = m. 

Problem 8. Check that the distribution f[3"91 is indeed an equilibrium solution of kinetic 
equation, i.e. I co11 = 0. 

If the reactions of annihilation of particle and antiparticle 



are in equilibrium, then from condition flHBl follows that /i 7 = and that the chemical 
potentials of particles and antiparticles are equal by magnitude and have opposite signs: 



If the equilibrium with respect to annihilation into two photons is maintained, while the 
annihilation into larger number of photons is out of equilibrium (such reactions are slower 
due to an extra power of the fine structure constant a and smaller phase space), non-zero 
chemical potential of photons can be developed. The observed CMB photons have zero 
or very small chemical potential \n/T\ < 10~ 4 . 

The equilibrium number density of bosons with = is: 



(39) 



b + b «-> 27 , 37 



(40) 



[i + \i = 0. 



(41) 



n b ^J 



(27T) 



MP) 




C(3)^T 3 /tt 2 w 0.12^T 3 , T>m; 
{2tx)^/ 2 g s {mT f I 2 e- m ' T , T < m , 



(42) 



where g s is the number of spin states. 

The number density of photons is equal to: 

n 7 = 0.2404T 3 = 412(T/2.728K) 3 cm 



-3 



(43) 
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where 2.728 K is the present day temperature of the cosmic microwave background radi- 
ation (CMBR). 

The equilibrium number density of non-degenerate (i.e. \x = 0) fermions is: 

_ / 3V4~0.09(7 S T 3 , T>m; 
f ~ \ {2Tx)^l 2 g s {mTfl 2 e~ m l T , T < m. [ ] 

The equilibrium energy density is given by: 

V — f dpp2E U*ft 



J 2tt 2 J exp[(£ - n)/T\ ± 1 

The total energy density of all species of relativistic matter with fi = is 

p rel = (tt 2 /30)^T 4 , (46) 

where g* = Yl [db + (7/8)gf] and gbj is the number of spin states of bosons or fermions. 
Problem 9. Calculate g* for T ~ 3 MeV. Answer: 10.75. 

Sometimes the total energy density is described by expression (146]) with the temper- 
ature depending number of species, g*{T) which includes contributions of all relativistic 
as well as non-relativistic species. 

The energy density of CMB photons is 

* = ^ a 02615 y^) 4 s ^ 4662 ■ 10-31 y^) 4 £5 • <47 > 

The relative contribution of CMB into the cosmological energy density is small but non- 
negligible: 

£W = 4.7-1(T 5 . (48) 

Heavy particles, i.e. those with m > T, have exponentially small number and energy 
densities if they are in equilibrium: 

If the annihilation stopped, the density of massive particles could strongly exceed the 
equilibrium one. If the particles are unstable with a large life-time then at later stage 
their distribution would return to the equilibrium one. 

Approach to equilibrium and deviations from it in homogeneous cosmology are de- 
scribed by the kinetic equation in FRW space-time: 

§ = (ft + Pd pi )fi = {d t - Hpta^fi = If 1 , (50) 

where p = —Hp and j^ co11 ^ is the collision integral, see below eq. f)62p . 

Problem 10. Why in the distribution function p and t are taken as independent 
variables, while above we treated momentum as a function of time, p = p(t)7 
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In terms of dimensionless variables: 



x = m a and jjj = Pjd 
the l.h.s. of kinetic equation takes a very simple form: 



Hx 



TCOU 



dx 



(51) 



(52) 



If the universe is dominated by relativistic matter, the temperature drops as T ~ 1/a 
and the Hubble parameter is expressed through X dbS'. 



H = 5.44 



9* mo 



10.75 x 2 mpi 



(53) 



In thermal equilibrium: g* = 2 for photons, g* = 7/2 for e ± -pairs, and g* = 7/8 for one 
family of left-handed neutrino. Since H = l/2t the relation between cosmological time 
and temperature of the primeval plasma has the form: t/sec ~ (MeV/T) 2 . 
For non-interacting particles, i.e. for I coU = 0, equation 



is solved as 



Hx^ = 
ox 



f = f(y,x in ) 



(54) 



(55) 



where initial value of the scale factor. Thus the distribution function main- 

tains its initial form in terms of variables x and y. For massless particles with non-zero 
chemical potential: / = feq{y,0i if they ever were in equilibrium. Here T ~ 1/a and 
£ = fi/T = const. Initially equilibrium distribution of massless particles maintains its 
equilibrium form even after interactions are switched off, as is observed in CMB. 

Let us check this important statement. The l.h.s. of the kinetic equation can be 
written as: 



(d t - Hpd p ) f, 



eg 







f E- fi 




Hp' 


dfeq 


[ T(t) - 




T T 


T 




dx 



(56) 



The factor in square brackets vanishes if (i — T/T = —H, which is true in the expanding 
universe, and if E(T/T) = —Hp which can be and is true only for E = p i.e. for m — 0. 

If the particle mass is non-zero and if the interaction is switched off at T ^> m, the 
distribution looks as an equilibrium one but in terms of p/T but not E/T. 

Problem 11. Find the distribution of massive particles decoupled at T < m. What if 
decoupling is non-instantaneous? 

In the equilibrium state and for vanishing chemical potentials entropy, S, in comoving 
volume is conserved: 



dS_ 
dt 



d 

It 1 



P 



P 



T 







(57) 



15 



In fact this equality is more general. It is true for any distribution function / = f(E/T), 
satisfying the condition of the covariant energy conservation, p = —3H (p + P) with ar- 
bitrary T(t). So we find: 



dt V T 



T \ T ) T 



(58) 



where the pressure is: 



(2tt) 3 3E j \T 



From this expression we can find P (remember that only T depends upon time) and 
integrate by parts to obtain: 

P=^(p + P), (60) 

which leads to the conservation law ( 1571) . 

Let us return now to the kinetic equation and specify the collision integral for an 
arbitrary process: i + Y <-> Z: 

T" = ^E/ dv z dvY5\p i + py-p z )[\A(Z^i + Y)\*l[fl[(l±f)- (61) 

4 Z,Y J Z i+Y 

|ii(t + y-,z)| 2 / i n/II( 1± /)]' 

Y Z 

where Y and Z are arbitrary, generally multi-particle states, Y\y f ^ s ^ ne product of phase 
space densities of particles forming the state Y, and 

** = Ilpfb < 62 > 

The signs '+' or '— ' in FJ(1 ± /) are chosen for bosons and fermions respectively. 
Equilibrium distributions by definition annihilate the collision integral: 

j«" [/(«*)] = (63) 

The standard Bose/Fermi distributions do that. It is easy to check that this is indeed 
true in T-invariant theory where the detailed balance condition holds, 

\A f (p)\ 2 = \A fl (p>)\ 2 (64) 

where p' is time reversed momentum, i.e. with the opposite sign of the space coordinate 
with respect to p. 

So after an evident change of variables in the collision integral we can factor out |^i/| 2 
and the integrand would be proportional to: 

n/ ire n(i ± f fm ) - n/ /in n(i ± f in ) = o (65) 
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It is easy to check that the functions f eq (1371) annihilate the collision integrals due to 
conservation of energy 

^E in = ^E fin (66) 

and if chemical potentials satisfy: 

This condition is enforced by reactions. 

Since we know that CP-invariance is broken and (mostly) believe that CPT invari- 
ance holds, we must conclude that the invariance with respect to time reversal, T- 
transformation, is broken as well. It means, in particular, that the detailed balance 
condition is invalid, |^4j/| 2 ^ \Afi\ 2 . Now a natural question arises: would the usual equi- 
librium distributions survive in T-violating theory? Let us check what happens with the 
collision integral for / = f eg . Due to eq. ( |65l) the integrand is proportional to 

hoii ~ n/ in (l ± f fm ) (\A lf \ 2 - \A fi \ 2 ) . (68) 

The last factor is non-vanishing if T-invariance is broken. However, due to S-matrix 
unitarity (or hermicity of the Hamiltonian) breaking of T-invariance is observable only if 
several processes participate and though each separate term is non-zero, the sum over all 
relevant processes vanishes [7J. 

It can be proven using S-natrix unitarity condition, SS^ = 1. If as usually we introduce 
scattering matrix: S — (I + iT), then it satisfies: 

i(T if - 4) = ~J2 T in T ] nf =J2 T l T -f ( 69 ) 

n n 

Summation over n includes integration over phase space. Instead of detailed balance a 
new condition of cyclic balance [7j 

[ dr k (\A kl \ 2 -\A lk \ 2 ) =0 (70) 

k 

ensures vanishing of I co u on / = f eq . Here dr k includes Bose/Fermi enhancement/suppression 
factors. For validity of this relation full unitarity is not necessary. Normalisation of proba- 
bility ^2fWif = 1 plus CPT invariance are sufficient. If CPT is broken, then the additional 
condition Ylf w if = 1> Ylf w fi = 1 would save the standard equilibrium statistics. How- 
ever, in the case that nothing above is true, equilibrium distributions would differ from 
the canonical ones. If the theory does not respects sacred principles of unitarity, hermic- 
ity, etc., the Pandora box of disasters would be open and equilibrium might deviate very 
much from the standard case or even not exist. 

6 Freezing of species 

As we discussed in the previous section, primordial plasma is typically in thermal equi- 
librium state in the early universe. When T dropped down to the so called decoupling 
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temperature, T^, (sometimes it is called freezing temperature, Tf) the interaction with 
plasma effectively switched off and the particles started to behave as free, non-interacting 
ones. There are two types of decoupling: 

1. Relativistic freezing, when Td > m. This is realised e.g. for neutrinos, or some other 
hypothetical weeakly interacting particles. Such particles by definition make hot dark 
matter (HDM) or warm dark matter (WDM). 

2. Non-relativistic freezing, when Td < m. Such particle make cold dark matter (CDM). 



6.1 Relativistic freezing 

Neutrinos decoupled at temperatures much larger than their mass, T ^> m v . It can be 
seen from the decoupling condition that the weak interaction rate became smaller than 
the expansion rate aw n < H at: 

G 2 F E 2 T 3 ~T 2 /Mpi, (71) 

where Gp = 1.166 x 10~ 5 GeV~ 2 is the Fermi coupling constant. 

Since the energy of relativistic particles is equal by an order of magnitude to the 
plasma temperature E ~ T we obtain: 

T f ~ m N (lO w m N /M Pl ) 1/3 ~ MeV. (72) 

More accurate calculations are necessary to establish if Tf is larger or smaller than m e , 
which is important for the calculations of the number density, n u , of relic neutrinos at the 
present time and for the cosmo logical bound on their mass, m„. 

For more accurate calculations of the decoupling temperature from the electron- 
positron plasma we will use the kinetic equation in Boltzmann approximation with only 
direct reactions with electrons, i.e. ve elastic scattering and z/P-annihilation taken into 
account: 

Ex = 80G F(9l + 9 2 R )y ^ (73) 

f u dx 37T 3 X 5 

where we define x = MeV/T. 

It is clear from this equation that the freezing temperature, Tf, depends upon the neu- 
trino momentum y = p/T, and this can distort the spectrum of the decoupled neutrinos, 
as we see in what follows. 

For the average value of the neutrino momentum, y = 3, the temperature of decoupling 
of neutrinos from is 

T„ e = 1.87 MeV, and T v ^ Vt = 3.12 MeV. (74) 

Problem 12. Derive equation ( 1731 . Find decoupling temperature of annihilation, vv^> e + e~ 
which changes the number density of neutrinos. 
Answer: T u& w 3 MeV and T v ^ Vt w 5 MeV. 

To take into account all reactions experienced by neutrinos including elastic ve and 
all vv scattering we need to make the substitution: (g\ + gp) — > (1 + g\ + gp) in eq. (175]) 
and find that neutrinos started to propagate freely in the universe when the temperature 
dropped below 

T Ue = 1.34 MeV and T v „ T = 1.5 MeV. (75) 
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6.2 Gershtein-Zeldovich(GZ) bound 

Consideration of thermal equilibrium and entropy conservation permitted Gerstein and 
Zeldovich to derive famous cosmological bound on neutrino mass [8]. Sometimes this 
bound is called Cowsic-McLelland bound but this is not just because paper [9] has been 
published 6 years after Gerstein and Zeldovich and contained a couple of inaccurate 
statements which resulted in overestimation of the bound by factor 22/3. 

We need to calculate the ratio n v j n 7 at the present time. The known from observations 
number density of photons in CMB, see eq. ( j43l) . allows to determine the cosmological 
number density of unobservable neutrinos. At neutrino decoupling the ratio is determined 
by thermal equilibrium: 

n v = n 9 = (3/8)n 7 (76) 

After decoupling n u was conserved in the comoving volume, i.e. n v a 3 = const but n 7 a 3 
rises due to e + e~-annihilation into photons. At first sight the rise of n 7 is difficult to 
calculate but entropy conservation ( 1571) makes the calculations trivial. 

After neutrino decoupling the number density of neutrinos fall down as cosmologi- 
cal volume, n„ ~ 1/a 3 . Before e + e~-annihilation but after neutrino decoupling, say at 
T ~ lMeV, the entropy of photons and electron-positron pairs was: 

S m ~(2 + 7/2)7f n a 3 n . (77) 

After annihilation it became: 

Sfin ~ 2T 3 in a 3 fin . (78) 

Since n 7 ~ T 3 , and Si n = Sfi n , the ratio of number densities of neutrino and photons 
dropped by the factor 4/11. If in the course of subsequent cosmological evolution the 
numbers of photons and neutrinos conserved in the comoving volume, the number density 
of neutrinos today must be 

3 

n v + U y = —n 1 = 112/cm 3 (79) 

This result is obtained under assumption of vanishingly small chemical potentials of neu- 
trinos, in other words for n v = n v . 

Energy density of neutrinos today should be smaller than the total energy density of 
matter, p m . This leads to the following upper bound on the sum of masses of all neutrino 
species: 

rn Vj < 94 eV Vth 2 . (80) 

Since h 2 w 0.5, Q m ~ 0.25, and masses of different neutrinos are nearly equal, as follows 
from the data on neutrino oscillations, we find m v < 5 eV. 

This limit may be further strengthen, if one takes into account that cosmological 
structure formation would be inhibited at small scales if SIhdm > 0.3 &cdm- Hence 
m u < 1.7 eV. Recent combined analysis of CMB and LSS leads to the bound: 

m^<0.3eV. (81) 
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For more detail and reviews see ref. [10] If the masses of neutrinos are close to this upper 
limit, their contribution into cosmo logical energy density at the present time would be 
non-negligible, Q u ~ 0.02, comparable to that of baryons, fi& ~ 0.04. 

One may argue that though massless neutrinos are 100% left-handed, i.e. they have 
only one helicity state, massive neutrinos have both spin states and hence their cosmolog- 
ical number density should be twice larger than calculated. However, it is not so because 
right-handed states did not reach equilibrium and their contribution may be neglected. 

Next question is how robust is the GZ bound. Is it possible to modify the standard 
picture to avoid or weaken it. The bound is based on the following assumptions: 

1. Thermal equilibrium between u, e ± , 7 at T ~ MeV. If the universe never was 
at T > MeV, neutrinos might be under-abundant and the bound would be much 
weaker. However, successful description of light element production at BBN makes 
it difficult or impossible to eliminate equilibrium neutrinos at the MeV phase in the 
universe evolution. 

2. Negligible lepton asymmetry. Non-zero lepton asymmetry would result in larger 
number /energy density of neutrinos plus antineutrinos and the bound would be 
stronger. 

3. No extra production of CMB photons after neutrino decoupling. Strictly speaking 
this is not excluded but strongly constrained. If the extra photons were created 
before BBN terminated, they might distort abundances of light elements. Late 
time creation of extra photons, after BBN, would lead to distortion of the energy 
spectrum of CMB and there is only very small freedom, not sufficient to change 
n u /n^ essentially. 

4. Neutrino stability on the cosmological scale, r v > tjj. If neutrino decays into another 
normal neutrino, e.g. /i M — > v e + X, the total number of neutrinos does not change 
and the limit on the mass of the lighest neutrino remains undisturbed, but heavier 
neutrinos are allowed. If the decay goes into a new lighter fermions, e.g. sterile 
neutrino, the bound may be weakened for all neutrino species. 

5. No late-time annihilation of v + v into a pair of (pseudo)goldstone bosons, e.g. 
majorons. For noticeable annihilation too strong coupling of neutrinos to majorons 
is necessary which is probably excluded by astrophysics. 

6.3 Distortion of neutrino spectrum 

As we mentioned above, see eqs. (I54H56I) . massless particles keep their equilibrium spec- 
trum even after the interaction is switched off. E.g. the spectrum of CMB photons is 
the equilibrium one with the precision better than 10~ 4 . However, it happened not to 
be true for neutrinos. The point is that neutrino decoupling is not instantaneous and for 
some time there coexist two components of plasma with different temperatures, weakly 
interacting with each other. Indeed, due to e + e~ annihilation the photon temperatures 
rises with respect to the neutrino temperature as, T 7 /T v « 1.4. 
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Due to residual interaction of neutrinos with hotter electrons and positrons some 
energy is transferred to colder neutrino sector. More energetic neutrinos decoupled later. 
As a result the spectrum is distorted [TTj : 

*/*//eff « 3 x 10" 4 - (^ - 3) (82) 

This analytical estimate was confirmed by precise numerical solution of the integro- 
differential kinetic equation; for discussion and the list of references see review [T2] . 
This effect leads to an increase of the effective number of neutrino species: 

AN V = 0.03 + 0.01. (83) 

The last 0.01 comes from plasma corrections [13], which diminish n e and n 7 with respect 
to unperturbed quantities at the same temperature. An increase of the number of neutrino 
species has negligible effect on BBN, ~ 10~ 4 , but may be noticeable in CMB measurements 
by the recently launched Planck mission. 

Note in conclusion of this section that though we mentioned above that neutrino 
temperature is approximately 1.4 times smaller than the temperature of photons, i.e. 
today it should be 1.95 K, would neutrino be massless, the distribution of neutrinos has 
the non-equilibrium form: 

/^[expfr/T^ + l]- 1 , (84) 

i.e the magnitude of neutrino momentum enters instead of energy and so the parameter 
T v does not have meaning of temperature. The correction f)82p is neglected here. 



6.4 Non-relativistic freezing 

If particles have sufficiently strong interactions, they decouple from primordial plasma at 
temperatures much smaller than their mass, Tf < m^. After that their number density 
stopped falling down according the the Boltzmann suppression law but remains constant 
in the comoving volume. The number density of heavy particles at decoupling is given by 

n h /n, « {m h /T f )V 2 e- m */ T f « 1, (85) 

so such particles may have masses much larger than permitted by GZ bound and can make 
cosmological interesting cold dark matter. The frozen number density of such particles is 
determined by the cross-section of their annihilation and is given by a simple expression, 
see e.g. [14] : 

^ {m \ m , (86) 

where m h /T f « ln((a ann v)m Pl m h ) ~ (10 - 50). 

We derive this expression in what follows. One can solve kinetic equation governing 
evolution of the number density of heavy particles numerically but it is instructive to make 
analytic calculations. Moreover, the results are pretty accurate. Analytic calculations of 
frozen abundances are usually done under the following assumptions: 
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1. Boltzmann statistics. It is usually a good approximation for heavy particles at 
T < m. 



2. It is assumed that heavy particles are in kinetic, but not chemical, equilibrium, i.e. 
their distribution function has the form: 

h = e~ E / T+ M , (87) 

where £ is the effective chemical potential normalised to temperature, £ = fi/T. 
Chemical equilibrium is enforced by annihilation which needs a partner whose 
number density is exponentially suppressed, while kinetic equilibrium demands en- 
counter with abundant massless particles. That's why chemical equilibrium stopped 
to be maintained much earlier than the kinetic one. 

3. The products of annihilation are in complete thermal equilibrium. 

4. Charge asymmetry of heavy particles is negligible and thus the chemical potentials 
for particles and antiparticles are equal, £ = +£. Annihilation would be much more 
efficient in the case of non-zero charge asymmetry. 

Kinetic equation under this assumption becomes an ordinary differential equation, which 
has been derived in 1965 by Zeldovich [15] and used for the calculation of the frozen 
number density of non-confined massive quarks in ref. [16]. In 1978 the equation was 
applied to the calculations of the frozen number densities of stable heavy leptons in 
ref. P~T] [T5] and after that it got the name Lee- Weinberg equation, though it would be 
more proper to call it Zeldovich equation. 

The equation has the following simple form: 

n h + 3Hn h = {a ann v){n h eq)2 - n 2 h ) , (88) 

where nh is the number density of heavy particles, n h eq ^ is its equilibrium value, and (a ann v) 
is thermally averaged annihilation cross-section multiplied by velocity of the annihilating 
particles: 

(crannv) = |^ J dp~ h dp~- h J dp ,dp / , 5\P in - P fin ) | A mn \ 2 e-( E * +E *>V T , (89) 

where dp = d 3 p/[2E (2vr) 3 ]. 

The integration in eq. ( 189]) can be taken down to one variable and we find j!9j : 



^ annV ^ = amSR^iV! / ds{s- AmDaannis^K! ) (90) 

° m h JX 2\ x ) J Am 2 h ^ h ' 

where x = rrih/T and s = (p + p) 2 . Usually x ^> 1 and a ann v — > const near threshold, so 
thermally averaged (av) is reduced just to the threshold value of av. The expression 
above can be useful if cross-section noticeably changes near threshold, e.g. in the case of 
resonance annihilation. 
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For derivation of eq. (1881) we start with the general kinetic equation: 



d t f-Hpd p f = I el + I ann , (91) 

where we take into account only two-body processes with heavy particles, l e \ and I ann 
are respectively collision integrals for elastic scattering and annihilation. At T < 
the former is much larger than the latter because of exponential suppression of the num- 
ber density of heavy particles, fh ~ exp(— m h /T). Since I d is large, it enforces kinetic 
equilibrium, i.e. canonical distribution over energy: 

f h = ex V [-E/T + m- (92) 

With such a form of fh we can integrate both sides of eq. (1911) over dp and the large 
elastic collision integral disappears, but the trace of it remains in the distribution f l92l . 
As the last step we express £(t) through rih- 

exp(£) = n h /n eq (93) 

and arrive to eq. (1HH1) . 

The equation can be solved analytically, approximately but quite accurately. Usually 
at high temperatures, T > m^, the annihilation rate is high: 

a ann n h /H > 1 (94) 

and thus the equilibrium with respect to annihilation is maintained, rth = n eq + 5n, where 
5n is small. It is convenient to introduce dimensionless ratio of number density to entropy 
r = nh/S, so the effects of expansion disappear from the equation: 

n + 3Hn = Sr . (95) 

Recall that S is conserved in comoving volume, eq. floT|) . 

By assumption r weakly deviated from equilibrium, so we can write r = r eq (l + 5r), 
where 8r <C 1- In this limit the solution of eq. ( l88l) can be found in stationary point 
approximation: 

Hxr' 

Since r eq exponentially drops down, 5r rises and at some moment 5r would reach unity. 
After that we will use another approximation, neglecting r^ q , integrate equation for r and 
obtain the final result ( 1861) . 

Another way to solve equation ( 1881) is to transform this Ricatti type equation to the 
second order Schroedinger type one and to integrate the latter in quasi-classical approxi- 
mation. 

Problem 13. Derive all above, in particular kinetic equation and solve it. 

Problem 14- Find frozen number densities of protons and electrons in charge symmetric 

universe. Answer: n p jn^ w 1CT 19 , n e jn^ m 10~ 16 . 

Problem 15. What number density would have anti-protons if (n p — n^)jn^ = 1CT 9 . 
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Let us apply the obtained results for calculation of the frozen number density of 
lightest supersymmetric particle (LSP) which must be stable if R-parity is conserved and 
is a popular candidate for dark matter. The annihilation cross-section is estimated as: 

uv ~ a 2 /m 2 s (97) 

Correspondingly the energy density of LSP would be: 

n 7 m| \n(a 2 Mpi/m s ) 
Psusy = m s n s « -* — s (98) 

Mpi 

For rris = 100 GeV, which is a reasonable value for minimal supersymmetric model, we 
find: 

Qsusy^ 0.05. (99) 

It is very close to the observed 0.25 and makes LSP a natural candidate for DM. 

Another interesting example is the frozen number density of magnetic monopoles, 
which may exist in spontaneously broken gauge theories containing 0(3) subgroup [2D] . 
The cross-section of the monopole-antimonopole annihilation can be estimated as 

(TannV ~ f/M^, (100) 

where Mm is the monopole mass. Correspondingly the present day energy density of 
magnetic monopoles would be [2*T] : 

a- M'j, . . 

Pm - 2n/r ■ (101) 

g 2 M Pl 

In fact slow diffusion of monopoles in cosmic plasma would slightly diminish the result 
but not much. 

If Mm ~ 10 17 GeV, as predicted by grand unified theories, the monopoles would over- 
close the universe by about 24 orders of magnitude assuming that their initial abundance 
was close to the thermal equilibrium one. This problem played a driving role for the 
suggestion of inflationary cosmology. 



7 Big bang nucleosynthesis 

BBN is one of the pillars of the standard cosmological model. It describes creation of light 
elements, 2 H, 3 He A He, and 7 Li, in the early universe when she was between 1 sec to 
200 sec old and the temperature ran in the interval from 1 MeV down to 60-70 keV. The 
calculated abundances of light elements are in a good agreement with the observation. 
This proves that our understanding of the universe when it was so young, is basically 
correct. 

The first stage of BBN is the freezing of the neutron-to-proton ratio, which determines 
the number density of neutrons for the second phase when formation of light elements 
took place. The n/p-freezing happened at T « 1 MeV and t ~ 1 s, while the light element 
formation occurred much later at T « 65 keV and t ~ 200 s. 
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The neutron-to-proton ratio is determined by the reactions: 

n + e + <-> p + v e (102) 
n + u e <-> p + e~ , (103) 

which frozen at T ^ 0.7 MeV (see below). After that r np remained almost constant, slowly 
decreasing due to the neutron decay: 

n <-» p + e~ + v e , (104) 

whose life-time is r n = 886 s. Since the formation of light elements started at t ~ 200 s, 
the decrease of r np due to decay was essential. However, the decay is not important for 
(n — p) freezing. 

It is more convenient to consider the ratio of neutron number density normalised to 
total baryon number density, r = n n /(n p + n n ), because latter is conserved in comoving 
volume at BBN epoch since baryonic number was conserved at low temperatures. 

Kinetic equation which governs the neutron to baryon ratio can be obtained from 
the general kinetic equation (!oT?|) with collision integral given by eq. (1621) in the limit of 
non-relativistic nucleons: 

. (1 + 3g 2 A )G 2 F 



2tt 3 



[A-(A + B)r] (105) 



where qa = —1.267 is the axial coupling constant of (n — p) weak current and the coeffi- 
cients A and B are 

i*oo poo 

A= dE 1/ Kf e (E e )[l-f u (E v )]\ Ee=Ei/+Am + dE e Kf p (E u )[l- f s (E e ]\ Ev=Ee+Am 

•JO J m e 

/■Am 

+ / dE e Kf p (E„)f e (E e ) \ E 

poo poo 

B = dE u K f u (E u ) [1 - f e (E e )} \ Ee=Ev+Am + / dE e Kf- e {E e ) [1 - f D {E u )} \ Ev=Ee+Am 

Jo J m e 

Mm 

+ / dE e K [1 - U(E V )\ [1 - f e (E e )} \ Ev+Ee=A m ■ (106) 
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where K = E^E e p e and we included terms describing neutron decay, the last ones in 
expressions for A and B. In precision calculations relativistic corrections and all form- 
factors of (n — ^-transformations are taken into account. 

The expressions for A and B take very simple form if e and v e are in thermal equilib- 
rium with equal temperatures. In this case A = B exp (— Am/T). Moreover, in essential 
range of temperature m e can be neglected and: 

B = 48T 5 + 24(Am)T 4 + 4(Am) 2 T 3 . (107) 

In this approximation equation (11051) can be easily solved numerically. 

We can make an estimate of the freezing temperature using the following simple con- 
siderations. The reaction rate versus Hubble rate is: 

T np (l + 3g 2 A )G 2 F B/27r 3 



H T 2 ^/0.6m Pl 
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The n/p freezing temperature can be approximately found from the condition T np /H = 1, 
that is 



and correspondingly (n/p)f = exp(—Am/T np ) « 0.135. Note that T np depends upon the 
number of neutrino species through , see discussion after eq. (I53p . 

When T drops down to T BBN — 60 — 70 keV, practically all neutrons quickly form 
4 He (about 25% by mass), 2 H (3 x 10~ 5 by number), 3 He (similar to 2 H), and 7 Li 
(10 -9 — 10~ 10 ). The calculated abundances span 9 orders of magnitude and well agree 
with the data. 

Question: why Tbbn is much smaller than nuclear binding energy, Ef, ~ MeV? The 
answer is below in this section. 

Almost all frozen neutrons, except for those which decayed before the onset of light 
element formation at T = Tbbn ~ 65 keV, form 4 He because of its largest binding energy 
equal to 7 MeV/nucleon. Correspondingly the mass fraction of 4 He can be estimated as: 



in a good agreement with the data. 

It is interesting that a small variation of the Fermi coupling constant would strongly 
change the amount of the produced 4 He and correspondingly the star properties. The 
stars might either have a deficit of helium or of hydrogen. 

Problem 16. Find the range of variation of Gp which is in agreement with (25 ± 1)% of 
the mass fraction of 4 He. Find the same for or the number of neutrino families. 

The light element formation proceeded through the chain of reactions: p (n, 7) d, 
d (^7) 3 He, d (d, n) 3 He, d (d,p) t, t (d, n) 4 He, etc. All reaction go through formation of 
deuterium, because due to low baryon density two body processes dominate. An absence 
of a stable nuclei with A = 5 results in suppression of heavier nuclei production. 

Let us turn now to the calculations of the temperature when the light elements were 
created. Naively one should expect this temperature to be close to the nucleus binding 
energy T ~ Ej,. However, a large number of photons make it possible to destroy the 
produced nuclei on the tail of their energy distribution, despite the Boltzmann suppres- 
sion. So one would expect that the temperature of nucleus formation is smaller than the 
binding energy by the logarithm of the ratio n^/n 7 . We will derive now the Saha equation 
and see that this is indeed that case. 

The equilibrium density of deuterium is determined by the equality of chemical po- 
tentials, jJ,d — pLp + fin- 
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Y = 2(n/p)/[l + (n/p)} w 24%. 
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where the chemical potential can be expressed through the number density as 




(112) 
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Thus in equilibrium we can express the number density of deuterium through the number 
densities of protons and neutrons: 



where Bp = m p + m n — vhd = 2.224 MeV is the binding energy of deuterium and n p = rjn^ 
with 7] = (3 = 6 • KT 10 < 1. Notation rj is used at consideration of BBN, while the same 
or similar quantity is denoted as f3 when baryogenesis is considered. 
The number density of deuterium, n d , becomes comparable to n n at 



where r/io = lO 10 ?]. At higher temperatures no (and densities of other nuclei) are much 
smaller than n n and thus Tbbn can be considered as the temperature when the formation 
of light nuclei started. Because of the exponential dependence on T the light nuclei 
formation proceeded during quite short time interval. 

Note, that due to the same effect, the hydrogen recombination took place at T ~ 0.1 
eV which is much smaller than the hydrogen binding energy, Eh = 14.6 eV. 

In Fig. [H taken from ref. [22], the calculated abundances of light elements as functions 
of the baryon-to-photon ratio are presented. 

Helium-4 slowly rises with rj, because for larger rj the moment of BBN becomes earlier, 
see eq. f )114p . and less neutrons decayed. 

Deuterium abundance quickly drops down with rising n because the probability of 
processing of deuterium to heavier more tightly bound nuclei, 4 ife, is larger with larger 
baryonic number density. So less deuterium survives at larger rj. High sensitivity of the 
primordial deuterium abundance to n allowed it to be the best way to measure the cos- 
mological amount of baryons before more accurate CMB measurements became available. 
This is why deuterium was called "baryometer" . 

Lithium-7 is formed in two competing processes with different dependence on 77. At low 
7710 < 3 the production predominantly goes through the reaction 3 if( 4 ife, 7) 7 Li. On the 
other hand, 7 Li can be destroyed by collisions with protons and with rising rj destruction 
becomes more efficient and 7 Li drops down. At larger r^ 10 > 3 the dominant process of 
creation is 3 He( 4 He,'y) 7 Be. Destruction of 7 Be by protons is less efficient because 7 Be 
has larger binding energy than 7 Li. So 7 Be production rises with rising rj. At lower T, 
when formation of atoms became non-negigible, 7 Be could capture electron and decay 
into 7 Li and neutrino. 

The BBN calculations are based on pretty well known low energy nuclear physics and 
theoretical uncertainties would not play a significant role in comparison of theory with 
observations, if we were able to observe these light elements at the epoch of their creation. 
However, we observe them now, while the results are obtained for very young universe, 
about 300 seconds old. So evolutionary effects should be taken into account. We will 
briefly describe the problems of comparison of theory with the data below. For more 
detailed discussion see reviews [2"2l |2"B"] . 

Helium-4 is very tightly bound nuclei and so it is not destroyed in the course of 
cosmological evolution. Hence the observed abundance of 4 He should be larger than the 
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Figure 1: The abundances of He 4 , D, He 3 and Li 7 predicted by the standard model of 
BBN. The bands show the 95% CL range. Boxes indicate the observed light element 
abundances (smaller boxes: ±2<r statistical errors; larger boxes: ±2cr statistical and 
systematic errors). The narrow vertical band indicates the CMB measure of the cosmic 
baryon density, while the wider band indicates the BBN concordance range (both at 95% 
CL). 
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primordial one. With the existing observation means 4 He can be observed only at low red- 
shifts in chemically evolved regions with the abundance which may be quite different from 
the primordial one. To deduce the primordial abundance of 4 He one needs to extrapolate 
to zero metallically. Namely the regions where A He is observed are contaminated by 
heavier elements. One can study the correlation of this elements with 4 He and extrapolate 
(linearly) the data to zero values of the metals (in astronomy all heavier than helium are 
called metals). Another source of uncertainty is not very well known fraction of ionised 
helium with respect to the total amount. This described by the so called ionisation 
correction, which is another source of uncertainty. 

Deuterium could be destroyed in the course of evolution in poorly controlled manner. 
Fortunately, in contrast to helium, the deuterium line can be observed at large red-shifts, 
z ~ 1, i.e. at the earlier stages of the cosmological evolution. If the clouds, where 
deuterium is observed, are not contaminated by heavier elements, there is a good chance 
that deuterium is primordial. However, the deuterium line is shifted from the hydrogen 
one only by 80 km/sec and the peculiar motion of the cloud may induce an essential 
systematic error. The average value of deuterium abundance is in good agreement with the 
value of 77 determined from CMB, but the individual values are rather strongly dispersed 
from 1.6 ■ 10~ 5 up to 3.5 • 10~ 5 . It would be interesting to understand the origin of such 
strong dispersion. 

Primordial lithium visibly creates a potential problem for BBN, but 7 Li is difficult to 
observe (in first generation stars) and maybe it is premature to worry about the disagree- 
ment. 

As a whole the data and theory are in a good agreement. Still there seems to be some 
"small clouds" . It would be very interesting if these clouds indicate new physics but most 
probably the resolution of the problems can be found in the traditional way when more 
accurate data are accumulated and better understood. 



BBN is sensitive to any form of energy which was present in the universe during light 
element formation. Indeed the universe cooling rate can be determined by equating two 
expressions for the cosmological energy density, namely the critical energy density and 
energy density of relativistic plasma with temperature T: 



compare to eq. ( 1531 ). The factor g* = 10.75 + 1.75 AN„ count the contributions from 
photons, e ± , neutrinos, and any other form of energy parametrized by AN V . Usually 
AN U is called the number of extra neutrinos, though the corresponding additional energy 
may have nothing to do with neutrinos. If AN U = 1, the additional energy is equal to the 
equilibrium energy density of one family of neutrinos plus antineutrinos with negligible 
mass (at BBN). In particular, if neutrinos are degenerate, i.e. their chemical potential \i 
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is non-zero, the additional energy density corresponds to 
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(116) 



where £ = fi/T. 

Problem 17. Derive eq. (11161) . Use eq. 045 1) and, if necessary, consult e.g. book [21], 
chapter V, sec. 58 for the calculation of the integral. 

There are two effects induced by variation of g*. First, larger AA^ leads to earlier n/p- 
freezing and higher n/p-ratio, see eq. (11091) . Second, with larger the BBN temperature 
(11141) would be reached faster and more neutrons could survive against decay prior the 
onset of the light element formation. Both effects work in the same direction and AA^ = 1 
would lead to an increase of A He by 5%. Depending upon the data analysis the existing 
observational limit is 



According to ref. [23], N u ~ 2.5 seems to be the best fit. What is it, a problem or an 
"experimental" error? 

Extra energy in electronic neutrinos have an additional and stronger effect on BBN 
because they can shift the equilibrium value of n/p-ratio: 



Hence the bounds on lepton asymmetries depend upon the neutrino flavour and are much 
more restrictive for i/ e , than for z/„ r : 



if a compensation between effects induced by and £ e is allowed. In absence of com- 
pensation the bounds are somewhat stronger. However, this results were obtained in the 
case of weak mixing between different neutrino flavours. In real case of large mixing angle 
solution to neutrino anomalies the bounds on all chemical potentials are equal and quite 
strong, see below, sec. M 

The results presented above are valid for the equilibrium distributions of neutrinos. If 
or v T were out of equilibrium at (n — p)-freezing their only effect would be a change 
of the energy density or, what is the same a contribution into AN U . As for u e , their 
deviations from equilibrium would directly change the freezing temperature, T np . The 
effect depends upon the distortion of the neutrino energy spectrum. In particular, an 
excess of v e at high energies would shift T np to lower values and would lead to a smaller 
n/p-ratio, which is opposite to the discussed above effect from an increase of the total 
energy density. 

Sensitivity of BBN to additional energy at T ~ 1 MeV is used for deriving bounds on 
the number density of new (light) particles or new interactions of, say, neutrinos. As is 
discussed above, the parameter g*, which describes the contribution of different species 
into cosmological energy density, cannot differ form the canonical value 10.75 more than 



AN U < 0.3 - 0.5. 
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by 1. If neutrinos are massive there should be additional right-handed states in addition to 
the usual left-handed ones. If right-handed neutrinos are created by the canonical weak 
interactions, their production probability is proportional to m 2 . The most favourable 
period for their production took place at T ~ 100 GeV through decays of W or Z bosons. 
The probability of production of vr is 

-n = ^ = 10 f^V T ^w + r z n z _ 



n u \ T J T 3 

One can check that vr have never been produced abundantly if m„ respects the GZ-bound. 

Right-handed neutrinos could be produced directly if there exist right-handed current 
induced e.g. by exchange of right-handed intermediate bosons, Wr. At some early stage 
of cosmological evolution vr might be abundantly created and this would endanger suc- 
cessful predictions of BBN. To avoid this problem vr should decouple before the QCD 
phase transition (p.t.). In this case their number and energy densities would be diluted 
by the entropy factor, as e.g. number density of the usual neutrinos were diluted by 
e + e~-annihilation discussed in sec. 16.21 The number of species above the QCD p.t. is 
58.25, which includes three quark families (u,d,s) with three colours and 8 gluons with 
two spin states plus the usual 10.75 from e ± , 7, and v. The energy dilution factor is 
(10.75/58.25) 4/3 « 0.105. That is if three v R had equilibrium energy density before QCD 
p.t., their energy density at BBN would make 0.3 of the energy density of the one normal 
neutrino. 

Since the ratio of the production rate of vr to the Hubble parameter is equal to: 

Tr/H = (T/Twf (m M /M WR f , (121) 
Vr would decouple before QCD p.t. if 

4/3 



m W R ( Tqqd \ 

m WL V 2 00MeVy 



(122) 



Here TV is the decoupling temperature of the normal neutrinos which is taken to be 3 
MeV. Thus mw R > (a few) TeV, which is an order of magnitude better than the direct 
experimental limit. 

When/if the precision in BBN will reach the level AN U < 0.1 — 0.2, the bound on 
mw R would be much stronger than above, namely about 10 4 TeV, because we will need 
to move to the electroweak phase transition or higher. 

Using similar arguments we can find a bound on the mixing between Wr and Wl' 

W x = cos9W L + sinO W R (123) 
The probability of production of vr at T < Tqcd in this case is given by: 

r = T R /H = sin 2 9 ( Tqcd ) " . (124) 
H/ V 2 00MeVy V ; 

The condition r < 0.3 results in sin 2 9 < 10~ 3 , or smaller. The effect does not vanish for 
mw 2 ~^ 00 • Due to the presence of v eR the {n — ^-transformation would be more efficient 
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and it would shift T np to smaller values, thus compensating the increase of g*, but the 
effect is small, ~ 2 . However, it may open window for very large 9 ~ 1. 

If neutrino has a non-zero magnetic moment, Vr would be produced in electromagnetic 
interactions: 

e ± + u L -> e ± + u R (125) 

Demanding AN V < 0.5 we obtain: 

\i v < 3 x 10~ 10 /iB (126) 

If there existed primordial magnetic field, then vr could be produced by the spin precession 
in this field and the bound on the magnetic moment would be: 

\l v < KTVb (Bprimord/G&USSy 1 (127) 

If there exist intergalactic magnetic fields with the strength Bi nt - ga i ~ 10 -6 Gauss and if 
these fields were generated in the early universe and evolved adiabatically, then: 

\i v < 10- 19 /iB • (128) 

If there exists mirror matter which is similar or identical to ours, BBN demands that 
the temperature of the mirror staff should be smaller than the temperature of the usual 
matter roughly by factor 2, see e.g. ref. [2"5] . 

More detail about the material of this and the next section can be found in review [12] . 



9 Neutrino oscillations in the early universe 

Neutrino oscillations in medium are modified in the same way as light propagation. It 
can be described by refraction index or what is the same (up to energy factor) by effective 
potential. In cosmological plasma the effective potential contains two terms |26j : 

V e a ff = ±C lV G F T 3 + C a 2 ^^ , (129) 

where Cj ~ 1 and rj is the plasma charge asymmetry: 

^ = 2 Vue + rj Vlt + t] Ut +r) e - n j2 (for u e ) (130) 
= 2^ + Vue + Vut - Vn /2 (for u„) (131) 

Effective potential is proportional to the amplitude of forward elastic scattering (the same 
as the optical refraction index) and is usually calculated in the lowest order in the coupling 
constant. In our case it is first order in Gf- 

In the local, 4-fermion, limit the Lagrangian describing elastic neutrino interaction 
has the form 

C v = (G F /V2)^g^ a (l + 7s )^ na (l+^ f , (132) 
/ 



32 



where / are fermions with which neutrino interacts (they include nucleons, charged lep- 
tons, and neutrinos) and is the proper coupling constant. This interaction form can 
be read-off any textbook on weak interaction. 

The first term comes from the averaging of the current J a '0z7a(l + 7s)^z over medium. 
Since the cosmological plasma is assumed to be homogeneous and isotropic, the average 
value of the space component of the current vanishes but (Jo) ^ 0, if the plasma is 
charge asymmetric. The result is proportional to the difference of the number densities of 
particles and antiparticles. This is the first term in eq. (11291) . The second term arises from 
non-locality of weak interactions due to W or Z boson exchange. At low energies the effect 
is proportional to q 2 /m^Y Z , where q is the momentum transfer. The second term looks 
formally as being of the second order on the coupling constant, but it is imply because 
the intermediate boson mass is written as m~ 2 ~ Gp/a, where a is the fine structure 
constant. In stellar interior V e ff is dominated by the first term, while in cosmology the 
second term takes over at T > 10 MeV. 

There is another significant complication in cosmology because in contrast to stars 
(except for supernovae at dense stage), the of neutrino absorption and production and of 
coherence breaking due to scattering are significant. Because of that the standard wave 
function description is not adequate, the system is essentially open and the density matrix 
formalism should be applied ■ The equation for density matrix has the following form: 



p=( -Hp—)p = i[H m + V eff ,p}+ dT(ujJ)(f l f I AA + --{p,ApA + }) + 



A and B are matrix amplitudes of scattering and annihilation respectively. The equation 
looks very complicated but numerical solution is possible. Moreover, in some cases, e.g. 
if there is MSW-resonance, even an accurate analytical solution can be found. 

Let us discuss now how neutrino oscillation may influence BBN bounds on lepton 
asymmetry. It seems that for large lepton asymmetry, L, the oscillations are expected to 
be inhibited due to large first term in effective potential (11291) . This term is non- negligible 
even if L ~ 10~ 9 and for L ~ 1 it might kill the oscillations at all. This is indeed the 
case for mixing between active and sterile neutrinos. However, for mixed active neutrinos 
off-diagonal terms in effective potential stimulate oscillations even in presence of large 



If the oscillations are not suppressed, individual leptonic numbers would not be con- 
served and e.g. initial asymmetry in, say, muonic sector would be equally redistributed 
between all three flavours: L e , L M , and L T . Thus the strict BBN bound on L e , eq. (11191) . 
became valid for all three leptonic numbers. The efficiency of redistribution of initially 
large lepton asymmetry by the oscillations was calculated in ref. [29], for more related 
references and analytical calculations see review [12J. It has been shown that for weak 
mixing (LOW solution for the solar neutrino anomaly) the redistribution of leptonic num- 
bers is insignificant, as is presented in fig. 2. However, for large mixing angle solution the 
equilibration of all leptonic numbers proceeded quite efficiently and they all became equal 
when the neutron freezing occured, see fig. 3. As a result, the following quite restrictive 
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bound on asymmetry of all neutrino flavours can be obtained [29] : 

16,1 < 0.07. (134) 

If this is the case, neutrino degeneracy cannot have a noticeable cosmo logical impact, 
in particularly on LSS and CMBR. The bound can be relaxed if neutrinos have new 
interaction with light Majorons. The potential induced by the majoron exchange inhibits 
early oscillations and different lepton numbers do not equalise [30] . 




Figure 2: Evolution of lepton charge asymmetry for LOW solution to solar anomaly. 



If there exists one or several sterile neutrino, z/ s , their mixing with the usual active 
ones would lead to several effects potentially observable in BBN: 

1. An excitation of additional neutrino species leading to positive AN„. 

2. Modification of spectrum of z/ e , because the oscillation probability depends upon neu- 
trino energy. 

3. Generation of large lepton asymmetry in v e . 

The impact of these processes on BBN was discussed in a large numbers of works, see 
review [T2j. In (non-realistic) case that active nus are not mixed and the mixing between 
v s and v a of a certain flavour is not resonance, the following bounds can be derived: 

(5m^ s /eV 2 )sin 4 2^ s = 3.16 ■ 1(T 5 ln 2 (l - AN V ) , (135) 
(Sml^JeV 2 ) sin 4 26 v » v ° = 1.74 ■ 10 -5 ln 2 (l — AN V ) . (136) 

They are noticeably stronger than the bounds obtained from direct experiment. 
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Figure 3: Evolution of lepton charge asymmetry for LMA solution to solar anomaly. 



If the realistic mixing between jy e ,fi,r are taken into account there are no simple analyt- 
ical expressions for the BBN bounds, but they have been obtained numerically in ref. [31] 
both in non- resonance case where they are similar to those above, eq. (113511136 1 ) . and for 
resonance case, where the bounds are much stronger. 



10 Inflation 

10.1 General features 

Inflation is a period of exponential (or more generally accelerated) expansion of the very 
early universe, with approximately constant Hubble parameter: 

a(t) ~ exp[Hit] . (137) 

It is the earliest time in the universe history about which we can say that it surely existed. 

Inflation is easy to realise e.g. by a scalar field, inflaton, with the energy-momentum 
tensor: 

T, u = d^<j> - (l/2)g, u [d a <j>d a <j> - U{ct>)} (138) 

If field slowly evolves, i.e. if (<9^0) 2 <C U (</>), then T^ u ~ g^ u and the vacuum-like equa- 
tion of state (127)) would be realised. 

The natural question, why do we need such a regime, has a very simple answer: 
inflation is the only known way to create the observed universe suitable for life. Inflation 
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naturally solves previously insolvable problems of Friedman cosmology: 

1. Flatness. Without inflation the fine-tuning of | O — 1 J should be 10~ 15 at BBN and 
10 -60 at the Planck era, see eq. (IHUj) Otherwise the universe would re-collapse in much 
less than 10 10 years, or expand too fast to make structures. Since \Q — 1| ~ exp(— 2Hiti), 
the duration of inflation should be: 

H i t i > 70. (139) 

More precisely the minimum duration of inflation may be somewhat smaller depending 
upon the heating temperature after inflation. 

2. Causality. The distance that CMB photons can propagate in the Friedman universe, 
before they stopped interacting, was calculated in sec. 13.21 eq. EH CMB photons stopped 
to interact with the surrounding medium after the hydrogen recombination at the red- 
shift z rec pa 10 3 , see eq. (11141) and below. At this moment the universe age was t pa 10 13 
s. Thus it is the maximal size of the region d c , which could be connected by interactions 
in Friedman cosmology. After recombination d c rises due to the cosmological expansion 
by z rec = 10 3 and today becomes ~ 10 16 s. The angular size of this path in the sky is: 

10 16 

e max = ^—~l°, (140) 

and regions outside this size should not know anything about each other. On the other 
hand, CMB comes practically the same from all the sky. At inflationary stage the causally 
connected region is exponentially large, so inflation could make all the observed universe 
causally connected if ifjtj > 70, the same as above. 

3. Inflation explains the origin of the initial push which induced cosmological expansion 
by antigravity of an almost constant scalar field, since for such a state P ~ —p and ac- 
cording the eq. (1T51) expansion speeds up and not slows down as naturally expected with 
attractive gravity. 

4. Inflation makes the universe almost homogeneous and isotropic at the present-day 
Hubble scale. Indeed, any perturbation with S/p/p ~ 1 with wave length A would trans- 
form into the perturbation with the same amplitude but with exponentially larger wave 
length: 

A -> exp(/fA)A (141) 

In other words, perturbations at fixed scale exponentially smooth down. 

5. Though inflation kills pre-existing perturbations, it creates its own small inhomo- 
geneities, at the level 10 -5 , but at astronomically large scales, which become seeds of 
large scale structure (LSS) formation. Inflation predicts adiabatic Gaussian density per- 
turbations with almost flat Harrison-Zeldovich spectrum. In terms of dimensionless grav- 
itational potential the spectrum does not contain any dimensional parameter and has 
approximately a simple form: <5\I/ ~ 5k /k. Deviations from flat spectrum agree with the 
data. 

Problem 18. According to the existing models, the natural duration of inflation is 
much larger than H{ti = 70, hence one should expect \Q — 1| negligibly small but in reality 
inflation predicts at horizon scale today: \Q — 1| ~ 10~ 4 . Why? 
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The idea of inflation was probably the most important breakthrough in cosmology 
of the XX century after the big bang one. Historically first paper where exponential 
expansion was invoked for solution of some problems of Friedman cosmology was that 
by Kazanas [32J, who suggested that exponential regime could solve the problem of the 
observed homogeneity and isotropy of the universe. A few months later a famous paper 
by Guth "Inflationary universe: A possible solution to the horizon and flatness problems" 
was published [33] • This work has initiated a stream of papers which remains unabated 
to the present day. In both scenarios the vacuum-like energy, which might dominate 
at first order phase transition, was suggested as a driving force of exponential expan- 
sion. It was soon understood that such mechanism was not satisfactory because it would 
create an inhomogeneous universe consisting of many relatively small bubbles in exponen- 
tially expanding vacuum-like background. The first workable mechanism of inflation was 
suggested by Linde [34J and Abrecht and Steinhardt [3H]. Probably the most beautiful 
inflationary mechanism, the so called chaotic inflation, was proposed by Linde [36J. For 
a review on these and discussed below issues, see refs. [37] . 

There was significant "pre-inflationary" literature directly related to the subject. The 
idea that the universe avoided singularity and underwent exponential period during which 
the mass of the cosmological matter rose by tens orders of magnitude was discussed by 
Gliner [3E] and Gliner and Dymnikova [39J . De Sitter like (exponentially expanding) non- 
singular cosmology was considered by Gurovich and Starobinsly [10] and by Starobin- 
sky jH]. In the last paper an important result was obtained that during "initial" expo- 
nentially expanding stage gravitational waves were produced which may be observable at 
the present time. If observed, it would be one of the strongest "experimental" indications 
to primordial inflation. 

Another prediction of inflation is the spectrum of primordial density perturbations 
which is already verified by the data. The pioneering calculations of the spectrum have 
been done by Mukhanov and Chibisov [12] and confirmed by many subsequent studies [37] . 

It was shown in the paper by Sato [I3J that exponential expansion induced by first 
order phase transition would never be terminated for certain under-critical values of the 
parameters. This happened to be a serious shortcoming of suggested later first inflation- 
ary scenarios. It was also notice by Sato [B] that exponential expansion might permit 
astronomically interesting antimatter domains. 

10.2 Models of inflation 

There are now many dozens of different mechanisms of inflation, see reviews [37j, and we 
do not have time to talk about them. We will discuss here only one which looks most 
economic and simple; it originated from the Linde's suggestion of chaotic inflation [36J. 

Let us assume that there existed in the very early universe a small piece where some 
scalar field (we call it inflaton) is almost constant. In other words all derivatives, d(p are 
small. We will see that this piece of the universe would expand exponentially and spatial 
derivatives would smooth further down. The size of this smooth part of the universe 
should be larger than the inverse Hubble parameter. It is difficult to evaluate probability 
of such a state but it may be unnecessary. If it existed (with whatever small probability), 
our large universe would evolve out of this microscopically small piece. 
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If spatial derivatives can be neglected, field satisfies the following equation of motion: 



+ 3#0 + U'(<f>) = (142) 

The second term takes into account cosmo logical expansion and is similar to the liquid 
friction in Newtonian mechanics. If H is in some sense large (we will specify below the 
proper conditions), then the equation can be reduced to the first order one: 

= -U'/3H. (143) 

Intuitively it is clear that in the case of large friction velocity is proportional to the force. 
This is the so called slow roll approximation which works pretty well in many inflationary 
scenarios. 

If the cosmological energy density, p, is dominated by slow varying inflaton field 0, 
then the Hubble parameter is equal to 



pi 



From this expression we can estimate the number of e-foldings while "lives" high in the 
potential U(4>) and the expansion is approximately exponential: 

N= [ H *=*fW&. (145) 



m? Pl J U'{<j>) 
For the power law potential, U(4>) = g<f) n , we find: 

47F "2 j,2 _ 47r 



N = r (0 1 - <f>j m » r n . (146) 

n m pl 1 n m pl 

For successful inflation N > 65 — 70 is necessary. It implies 0j n > 2.5n 1 ^ 2 m P i. At first sight 
it looks disturbing if something is larger than m v \. However this is not the case to worry 
about, because the observable quantity is the energy density of and it would remain 
much smaller than m Pl because of small mass of and the self-coupling constants, as 
we see in what follows. 

For the validity of the slow roll approximation the following two conditions are to be 
fulfilled: 

< 3#0. (147) 

and 

2 < 2f/(0) (148) 

These conditions are satisfied if 



If" , 87T 

77" « 2~ 149 
U 3m Pl 



38 



E.g. for massive free field U = m 2 2 /2 (harmonic potential) the slow roll approximation 
would be valid if: 

<p 2 > (4?r/3) m 2 Pl (150) 

With (ft exactly at the lower limit the number of e-foldings is not enough but a slightly 
larger would do the job. The harmonic potential would not exceed the Planck value if 

2 < m 4 Pl /ml (151) 

If we take equal to the upper bound, (p in = mpjm^, and ~ 10~ 6 m P i, which is 
demanded by the condition of sufficiently small density perturbations, the number of e- 
folding would be huge: iV = 10 13 . The characteristic time when all this happened is tiny, 
t inf ~ 1(T 31 s. 

Problem 19. Find N and t in j for U = A0 4 , assuming A = 10~ 12 and initial U (0) = m Pl . 



10.3 Particle production by inflaton 

During inflation the curvature of the inflaton potential should be smaller than the Hub- 
ble parameter, U"(<f)) < Hi, see the slow roll conditions above. In this regime inflaton 
monotonically but slowly moved down to the minimum of U{<p). Simultaneously decreases 
the Hubble parameter. At some moment the second derivative in eq. (I142p became non- 
negligible and started to oscillate near minimum with frequency ui > H. The expansion 
regime changed from the exponential to matter dominated one, if U(4>) = m 2 i 2 /2, or to 
relativistic regime, if U(<f>) = A0 4 /4. 

Problem 20. Solve eq. (11421) with potential U = m 2 <j) 2 and find cosmological scale 
factor a(t) and H(t) assuming that the cosmological energy is dominated by <fi. Use 
expression (j!38j) for the energy-momentum tensor of <ft. 

As is well known, a time varying field produces particles to which it is coupled and 
especially efficiently those whose mass is smaller than the frequency of the oscillations. It 
is exactly what happened with the inflaton. Empty, cold universe filled with oscillating 
4> exploded with creation of hot relativistic particles. It is almost as is described in the 
Bible: "Let there be light". This moment is proper to call big bang, though it may be 
not generally accepted terminology. 

Particle production by inflaton have been first calculated perturbatively in ref. [4"5] . 
It was shown that due to weakness of inflaton coupling to matter, perturbative particle 
creation is rather slow and the temperature of the universe heating after inflation is 
relatively small, much smaller than p 1 /^, where pfi n is the inflaton energy density to the 
end of inflation. Let us consider an example when the inflaton is coupled to fermions 
through the Yukawa coupling: 

Um = g<Kt)W (152) 
where <f)(t) is supposed to be classical field satisfying eq. ( 11421) : 

m = rn_ Pl sinm,(t + t ) ^ 
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In the lowest order of perturbation theory the amplitude of production of pair of fermions 
with momenta k\ and k 2 is 



A(h, k 2 ) = g [ d 4 xa(t)(k 1: k 2 \^\0) = 4==<K k i + k a )ii(fciMW(£i + E 2 ), (154) 
J V4-&1-&2 

where the standard decomposition of ij) and ip in creation-annihilation operators is used, 
u and v are the Dirac spinors, Ei — \ki\ is the particle energy, and 

4>{u) = J dte lut (j){t). (155) 

The probability of particle production per unit volume is 

N ^V=vJ W^ 1 - 41 = £ L ^ ^ (156) 

The volume factor V as usually comes from the square of 5-function: 

[5(kx + k 2 )] 2 = 5(k! + k 2 )y/(27r) 3 (157) 

If uj ^> t~ x then the integration over time in the interval At 3> uj" 1 in eq. (11551) gives 
5(u — 2E) and the square of it would be 5t5(u — 2E)/(2ir). So we obtain for the rate of 
the fermion production per unit time and unit volume: 



Nf m\uj 2 9 2 m 2 pi 
At ~ 8?r ~~ 24vr 2 (t + t ) 2 



N f = 4 L = ^±=- = - V. " (158) 



This corresponds to the following decay rate of the field 0: 

Nf Q 2 

r « = wrk m «- (159) 

This is, as one can expect, the decay width of the 0-meson. 

Because of the particle production cj){t) should decrease faster than it given by eq. fll53[) 
For this can be taken into account by the substitution <f>(t) — > <p{t) exp(— T^t). 

The rate of thermalisation of the produced fermions is as a rule larger than the expan- 
sion rate. In this case the temperature of the plasma can be evaluated as follows. Assume 
that the particles are produced instantly at the moment when the Hubble parameter 
H = l/2t becomes equal to the decay rate and that t t . The energy density of the 
produced fermions is 

Pf = = ( 16 °) 

and correspondingly the temperature of the universe heating is 



V4 , go X l/4 



\ir 2 gj V128tt 5 ^ 



g^/mpim^. (161) 



10 



For a more accurate evaluation of let us take into account non instant character 
of the particle production and the decrease of the amplitude of the oscillations <f>o(t) not 
only because of the Universe expansion but also due to the particle production. With 
these factors taken into account the energy density of the produced fermions satisfies the 
equation 

p s = T m -m Pf , (162) 

where Y $ is given by eq. (11591) and p$ = m^,0o(t) is the energy density of oscillating (f)(t). 
We assume also that the total energy density is equal to the critical one: 

"* + "< = h m " H2 = ^w (163) 

In the last equation it is assumed also that the non-relativistic expansion law is valid, 
that is > pf. Thus eq. (I162p takes the form: 

p f = T <i> p c -(4H + T (f> )pf. (164) 

Integrating it with the initial condition p/(0) = we obtain for the early MD-stage: 

F 2 2 p -rt r rt 
p = \ m Pi e / dxe x x V3 (i 6 5) 

Pf 6vr[r(t + t )]8/3 J 1 J 

This equation is valid till the energy density of the produced fermions becomes non- 
negligible, e.g. till pf ~ p^ « Pc/2. This is realised at Tt = 1.3 At that moment 
P f pa g A m 2 ^m 2 pl /hQ0ix i and the temperature is 

/ 3 \V 4 

Th ~[50^iJ 9 ^^ Fl (166) 

This is about twice smaller than result (I16ip . 

In many interesting cases perturbative approximation to particle production do not 
adequately describe physics of the process. Non-perturbative calculations have been pio- 
neered in ref. jl6] and further developed in ref. [17] with an emphasis on the parametric 
resonance enhacement. Due to that the inflaton decay can proceed much faster than ex- 
pected from perturbation theory and the initial particle state might be far from thermal 
equilibrium. In particular, heavy particles with mass larger than would-be temperature 
could be produced. It was also noticed [IB] that production of heavy particles by gravi- 
tational field at the final stage of inflation might be essential as well. 

We conclude here that the heating temperature after inflation is model dependent 
and most probably is not large, Th < Equt ~ 10 15 GeV, i.e. GUT era probably never 
existed. It solves the problem of overabundant magnetic monopoles. Initial hot universe 
might be far from thermal equilibrium and very heavy particles could be produced by the 
cosmological gravitational field |48j . 
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10.4 Inflationary density perturbations 

The increase of quantum fluctuations of scalar field in the exponentially expanding space- 
time gives rise to density inhomogeneities in the Universe. Physically this phenomenon is 
connected with different moments of the end of inflation in different space points (in an 
appropriate reference frame) due to small spatial fluctuations of the inflaton. Below we 
will follow the simple presentation of review [49J. For detailed rigourous treatment see 
book [SD] 

Exponential expansion transforms microscopically small wave lengths of quantum fluc- 
tuations into astronomically large ones and so produce natural candidates for initial den- 
sity inhomogeneities which could be the seeds of the large scale structure formation. In a 
sense this task is even over-fulfilled because the inhomogeneities 5p prove to be too large 
for natural values of the parameters of the theory. In particular to get Sp/p ~ 10~ 4 on 
the galactic scale the self-coupling of the inflaton should be A < 1CT 12 . 

Since we expect that the density fluctuations have originated from quantum fluctuation 
during inflation, we need to say a few words about scalar field quantisation in De Sitter 
space-time. We assume that the field is massless because we are interested in <C H. 

Field (f) is quantised in the De Sitter space-time along the same lines as in flat space- 
time. The starting point is the expansion in creation-annihilation operators: 

/d z k 
WF^ [ate ~' txm + h - c] ' (167) 



where ujk = Vk 2 + m 2 and 4>k{t) satisfies the Furrier transformed free equation of motion 
(in curved space-time): 

k 2 

<f) k + 3H<p k + —<p k = 0. (168) 
cr 

Operators a k and a^, obey the standard commutation relation 

[a kl a+} = 5 {3 \k-k'). (169) 

For massless scalar field the solution of eq. (11681) is expressed through the Bessel functions: 

Mt) = C lk {krfl 2 [J^ /2 {kr) + C 2k J 3/2 (kr)}, (170) 

where r = exp(—Ht)/H is the conformal time and the coefficients Ci k are determined 
by the matching <p k to the corresponding expression in flat space-time, <p k —>■ exp(—iu k t) 
for H — > and r — > H~ l — t. Substituting asymptotic values of the Bessel functions for 
kr — » oo into eq. (I170p we find 

\j 2 k 

C 2k = -i, (171) 

where a is a constant phase. 
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In the limit of large t (or r — > 0) one easily finds: 

\M 2 « (i + ^- 2i ") ("2) 

Let us extract from the inflaton field <f)(x,t) classical homogeneous part (po(t): 

(j>(x,t) = <j>o(t) + 6<i>(x,t), (173) 

where 5<p(x, t) describes small quantum fluctuations. On inflationary stage it satisfies the 
equation: 

<ty + 3H6j> - e~ 2m d 2 H - d2 ^[f o) 6<f> = 0. (174) 

o<p l 

For large Ht the third term in the equation can be neglected and S(p satisfies the same 
equation as <fio does, see eq. (11421) . Equation (11741) has two solutions. One of them 
decreases as exp(— 3Ht) for d 2 V/d<p 2 <C H 2 . The second solution varies relatively slowly. 
So at large t the first solution can be neglected and we can write 

6<f>(x,t) = -5r(x)j>o(t). (175) 

If is small this is equivalent to x-dependent retardation of the classical field motion to 
the equilibrium point: 

<f>(x,t) = <p (t - 5t(x)) . (176) 

Correspondingly inflation ends at different moments in different space points. This is 
the physical reason for generation of density perturbations. Since the energy density 
in the universe during inflation is dominated by the inflaton field <p, one can write 
p(x,t) = p(t — Sr(x)) forgetting possible subtleties connected with the freedom in co- 
ordinate choice. Thus we get 

^ = -5r^ = AH5t(x). (177) 

At the last step the relation p/p = —4H has been used. It is valid on RD-stage which by 
assumption was formed in the heated universe when inflation was over. It is convenient 
to make the Furrier transform with a different integration measure: 

^rlA<KM)| 2 . (178) 

The density fluctuations can be now expressed in terms of A0fc as 

6 -l = AH^. (179) 

P 00 

The Fourier amplitude of the fluctuations is evaluated with the help of eqs. ( 1170ll72j) as: 

^A = ^{^W e ' 2m T ■ (180) 
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Recall that this expression is valid for massless free field. Since S(p(x, t) satisfies eq. 
( 11741) the approximation is valid till k 2 exp(— 2Ht) > \d 2 U/d<p 2 \. On the other hand, the 
proportionality condition (I175p is valid in the opposite limit when one can neglect the 
decreasing part of the solution. So for the evaluation of the density fluctuations one has 



to substitute expression ( 11801) in the boundary region: 

One should remember that U{4>) depends on t through <p = <fio(t) so eq. (I18ip implicitly 
determines tb- Finely we find: 



Ap H 2 



P 7t3/20 o [ tb (k)) 



i d 2 v(Mh(k))) 

H 2 deb 2 



1/2 

(1821 



This expression for the fluctuation spectrum is valid till the fluctuation wave length reaches 
the horizon. This is indicated by sub-h in the left hand side. 

Since <f)(t) is a slowly varying function of time, the fluctuation spectrum weakly de- 
pends on k. This flat spectrum is known as the Harrison-Zel'dovich spectrum For a 
satisfactory description of the universe structure the fluctuations should have the value 
Ap/p\ h « lO" 4 . 

As an example let us consider density perturbations in the model with the potential 
U(4>) = — A0 4 /4. This potential is not bounded from below, so we can trust it only for 
sufficiently low <fi. Homogeneous classical field 4>o(t) satisfies the equation: 

+ 3H<p - A0 3 = 0. (183) 

We assume that the term can be neglected, see subsection 110. 2\ and check the validity 
of this approximation on the explicit solution. The latter has the form: 

1 /2 

W)=(^\ (*/-r V2 , ( 184 ) 



V 2A 

where tf is approximately the moment when inflation ended. In a realistic model U(4>) 
is bounded from below, so <fr(t) does not tend to infinity. Still solution (11841) gives a 
satisfactory approximation up to t slightly less than tf. After that moment the rise of <p 
should turn into oscillations around the equilibrium point. 
The neglect of <fi is justified if the condition 

< 1 (185) 



3#0 2H(t f - t) 
is fulfilled. Substituting solution ( 11841) into eq. (1182ft we arrive to 

1/2 



^l-(^) (186) 
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where in accordance with eq. A181 j) : 



k 1 2H(t f -t b ) . , 

Ht b = ]n- + - In l/ Q b> . (187) 

The second term in this expression is evidently small. Ht b can be approximately expressed 
through the comoving wave vector k and the corresponding to it present-day physical scale 
/o as: 

ka(tf) kT Q 

where a(t) is the scale factor, To is the present-day value of the CMB temperature, and 
Th is the heating temperature at the end of inflation. It is taken into account that 
the physical momentum changes inversely proportionally to the scale factor. In particular 
during inflationary stage p = k exp(-Ht) and during Friedman stage p basically decreases 
as inverse temperature. 



As a result 



and finally we obtain 



™, = lnf^ (189) 



( 8 \ ,1/2 i 3/2 ' 



p h \3tt 3 J b 



A^ln 3 / 2 ^, (190) 



where b = (2itTh/ HTq). For Th = 10 15 GeV which is possibly rather high, but not 
unreasonable, we obtain b ~ 15m ~ 1.5 • 10~ 15 years. Hence the density fluctuations on 
the galactic scale Zo = 10 6 years are 

Ap/p « 10 2 A 1/2 (191) 

So that A should be tiny, A m 10~ 12 , to give rise to a proper value of the fluctuations, 
Ap/p w 10~ 4 . This is a common shortcoming of inflationary models. Such a small 
value means in particular that the inflaton should be a gauge singlet or more complicated 
scenarios are necessary. Otherwise the interactions with gauge bosons would generate too 
big effective coupling, A e //0 4 with A e // ~ « 2 ~ 10~ 4 . 

The perturbation spectrum is not exactly flat but slightly deviates from the Harrison- 
Zeldovich one. This is a general feature of all inflationary models. 

Generation of gravitational waves during inflation can be considered along the same 
lines. Moreover, the equation of motion of spin eigenstates of gravitons coincides with 
the equation of motion of massless scalar field. It may be instructive to note that gravi- 
tational waves are not created in eternal De Sitter background. One can check that the 
Bogolyubov coefficients which describe particle production are trivial. But when expan- 
sion is changed from exponential to, say, power law, the production of particles and waves 
becomes possible. 
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10.5 Inflationary conclusion 

Inflation seems to be practically an experimental fact. It nicely explains the observed fea- 
tures of the universe. In particular, the observed spatial flatness today, D, — 1 ± (~ 10 -2 ) 
and very exact fine-tuning at, say, BBN, Q = 1 ± 10 -15 , naturally arises because of expo- 
nential expansion in the early universe. The observed spectrum of density fluctuations is 
predicted to be nearly flat, but not exactly flat. This is another successful prediction of 
inflation. The only "missing link" is not yet observed, namely, long gravitational waves, 
which may be accessible to LISA. However, one should keep in mind that though an ob- 
servation of such waves will be a strong argument in favour of inflation but if they are not 
observed, inflation still will not be killed, because the amplitude of inflationary gravita- 
tional waves is model dependent and may be so small that they will escape observation 
by the near future antennas. 

11 Baryogenesis 

The observed part of the universe is 100% C(CP)-asymmetric. Up to now no astronom- 
ically significant objects consisting antimatter have been detected. There is only matter 
and no antimatter, except for a small number of antiprotons and positrons most probably 
of secondary origin. From the bounds of the flux of 100 MeV gamma rays one can conclude 
that the nearest galaxy, if dominated by antimatter, should be at least at ~10 Mpc [51] . 
However we cannot say much about galaxies outside of our super-cluster. Observed col- 
liding galaxies at any distance or galaxies in the common cloud of intergalactic gas are 
of the same kind of matter (or antimatter?). In particular, the fraction of antimatter in 
two colliding galaxies in Bullet Cluster is bounded by n B /riB < 3 x 10~ 6 [52]. In charge 
symmetric universe the nearest antimatter domain should be practically at the cosmolog- 
ical horizon, Is > Gpc, because of very efficient annihilation at an early stage [53] • Still 
smaller clumps of antimatter are allowed in our neighbourhood. 

However, one should bear in mind that these bounds are true if antimatter makes 
exactly the same type objects as the observed matter. For example, compact objects 
made of antimatter may escape observations and be quite abundant and almost at hand. 

A natural question arises in this connection: is matter predominance accidental (a re- 
sult of asymmetric initial conditions) or dynamical? Inflation proves that it is dynamical, 
originated from some physical processes with non-conserved baryon number. Sufficient 
inflation should last at least (~ 70 Hubble times) with practically constant Hubble pa- 
rameter. This could happen only if the energy density was approximately constant, see 
eq. ({TBI . However, if baryons are conserved, the energy density associated with baryonic 
number) cannot be constant and inflation could last at most 4-5 Hubble times. Indeed if 
baryonic charge were conserved then it would remain constant in the comoving volume, 
i.e. B ~ 1/a 3 . The energy density of bearers of this baryonic number cannot stay con- 
stant as well but should evolve as ps ~ 1/°™; where n = 3 for non-relativistic matter and 
n = 4 for relativistic matter. 

Using the observed value of the cosmological density of the baryons and assuming its 
conservation, we find that at the early relativistic epoch (RD), e.g. at BBN the energy 
density related to this baryonic number should be p B ~ 10~ 7 p to4 . During evolution at 
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RD-stage the ratio Pb/ Ptot remained practically constant. Let us go backward in time 
till inflation. At inflation p tot m const, but baryons excluded. The energy density of the 
latter, if their number is conserved, should evolve as pb ~ exp(—4Ht). It means that 
4-5 Hubble times back sub-dominant baryons were dominant and ptot ~ Pb, which could 
not stay constant more that 4-5 Hubble times. This is surely insufficient for solution of 
cosmological problems discussed in the previous section. 

It was suggested by Sakharov in 1967 {S3] that the cosmological baryon asymmetry 
could be generated dynamically in the early universe if the following three conditions are 
fulfilled: 

1. Non- conservation of baryons 

2. Breaking of C andCP symmetries. 

3. Deviation from thermal equilibrium. 

Note in passing that none of them is obligatory - we present examples below. But first 
let us discuss normal baryogenesis and validity of these three conditions. 

I. Non-conservation of baryons is justified theoretically. Grand unified theories, non- 
minimal SUSY, and electoweak (EW) theory predict that baryonic number is not con- 
served, AB 0. However, at the present time this prediction is not confirmed by direct 
experiment. Despite an extensive search, only upper bounds on proton life-time and pe- 
riod of neutron-antineutron oscillations are established. The only "experimental piece of 
data" in favour of baryon non-conservation is our universe: we exist, ergo baryons are 
not conserved. Half of century ago from the same experimental fact, our existence, an 
opposite conclusion of baryon conservation was deduced. Theory is an important input 
in understanding of what we observe. 

II. C and CP violation are discovered and confirmed in direct experiments. At the first 
part of the XXth century the common belief was that physics was invariant with respect 
to separate action of all three transformations: mirror reflection, P, charge conjugation, 
C, and time reversal, T. The weakest link in this chain of discrete symmetries was P, 
found to be broken in 1956 [SB] . 

It was immediately assumed that the world was symmetric with respect to the com- 
bined transformation from particles to mirror reflected antiparticles, CP. Both P and C 
are 100% broken in weak interactions but still some symmetry between particles and an- 
tiparticles was saved. This symmetry crashed down pretty soon, in 1964 [56]. After this 
discovery life in the universe became possible. 

Why CP-breaking is necessary for generation of cosmological baryon asymmetry but 
C-breaking is not enough? A formal answer to this question is the following. Let us 
assume first that C is conserved and that the universe was initially in C eigenstate, i.e. 



where \u) is the wave function of the universe and \rj\ = 1 is a constant. This means, 
in particular, that the universe had initially all zero charges because charge operator 
anticommute with C-transformation. May some non-zero charge, e.g. B, be generated 
dynamically? The answer is negative because due to C-invariance the Hamiltonian of the 
system commutes with C-operator, [C, TC] = 0. The time evolution of B is given by: 



C\u) = T)\u) 



(192) 



B(t) = (u\e- mt j£e mt \u). 



(193) 
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Let us insert into this equation the unity operator I = C l C: 

B(t) = (u\Ie~ mt IJ*Ie mt I\u} = -B(t), (194) 

taken that C J®C~ X = — Jq . Thus in C-conserving theory B(t) = B in = 0. 

The same arguments with CP instead of C prove that charge asymmetry cannot be 
generated, if CP is conserved and the universe is an eigenstate of CP: 

CP\u)=7]\u). (195) 

In rotating universe charge asymmetry might be generated even if CP is conserved. 
Global rotation can be transformed into baryonic charge! However, it seems difficult 
to realise such an idea. New long-range interactions, possibly very unusual, are needed. 

For B-generation in elementary (local) processes no assumption about the universe 
state is necessary: if CP is conserved, no asymmetry is generated through particle decays 
or reactions. We will discuss this below in concrete examples. 

At the present time only CPT-symmetry survived destruction. It is the only one 
which has rigorous theoretical justification, CPT-theorem, based on solid ground: of 
Lorenz-invariance, canonical spin-statistics relation, and positive definite energy. Still 
models without CPT are considered, e.g. for explanation of some neutrino anomalies and 
for baryogenesis. 

III. Thermal equilibrium is always broken for massive particles, but usually very little. 
To estimate the effect let us approximate the collision integral in kinetic equation (I50p as 
Icon = r(/ e? — /), where T is the interaction rate. Let us assume that V is large so that 
the deviation from equilibrium is small, 6f/f eq <C 1. Substituting into the l.h.s. of kinetic 
equation f eq , as is done in eq. fl56l) . we find: 

5f Hm 2 Tm 2 , . 

— — — pa ^ (1961 

f eq TTE TEmpi 1 } 

Since the Planck mass is very large the deviation from equilibrium might be significant 
only at large temperatures or tiny T. However, if the fundamental gravity scale is near 
TeV [57], equilibrium could be strongly broken even at the electroweak scale. 

Another source of deviation from equilibrium in the cosmological plasma could be 
first order phase transition from, say, unbroken to broken symmetry phase in non-abelian 
gauge theories with spontaneous symmetry violation. There might be a rather long non- 
equilibrium period of coexisting two phases. 

There are plenty scenarios of baryogenesis each of them one way or other performing 
a rather modest task explaining only one number, the observed asymmetry: 

= (n B -n B )/ nj = 6- 1(T 10 , (197) 

found from the analysis of two independent measurements: of light element abundances 
created at BBN and of angular fluctuations of CMB. 

It is a great challenge to astronomers to check if (3 is constant or it may vary at different 
space points, (3 = (3{x). What is characteristic scale Ib of variation of baryonic number 
density? May there be astronomically large domains of antimatter nearby or only very far 
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away? Answers to these questions depend, in particular, upon mechanism of CP violation 
realized in cosmology, which are described below. For more detail see lectures [5H]. There 
are three possibilities for CP-breaking in cosmology: 

1. Explicit, realized by complex coupling constants in Lagrangian, in particular, com- 
plex Yukawa couplings transformed by the vacuum expectation value of the Higgs field 
((p) ^ into non-vanishing phase in CKM-mixing matrix. However, in the minimal stan- 
dard model (MSM) based on SU(3) x SU(2) x 17(1) CP-violation at T ~ TeV is too 
weak, at least by 10 orders of magnitude, to allow for generation of the observed baryon 
asymmetry. Indeed, CP-violation in MSM is absent for two quark families because the 
phase in quark mass matrix can be rotated away. So at least three families are necessary. 
It could be an anthropic explanation why we need three generations. 

If masses of different up or down quarks are equal, CP violation can be also rotated 
away because the unit matrix is invariant with respect to unitary transformations. If the 
mass matrix is diagonal in the same representation as flavour matrix, CP-violation can 
also be rotated away. Thus CP-breaking is proportional to the product of the mixing 
angles and to the mass differences of all down and all up quarks: 

A_ ~ sin #i2 sin #23 sin #31 sin 5 (m 2 — m^)(m t 2 — m 2 ){m 2 — ml) (198) 

(ml - m 2 s )(m 2 b - m 2 d )(m 2 s - m 2 d )/ M 12 . 

At high T > TeV, where electroweak baryon nonconservation is operative, the character- 
istic mass M ~ 100 GeV and A_ ~ 10~ 19 . So for successful baryogenesis an extension of 
MSM is necessary. 

2. Spontaneous CP violation [59], which could be realized by a complex scalar field <3> 
with CP-symmetric potential, with two separated minima at ($) = ±/. The Lagrangian 
is supposed to be CP-invariant but these two vacuum states have the opposite signs of 
CP-violation. Such CP-breaking is locally indistinguishable from the explicit one but 
globally leads to charge symmetric universe with equal amount of matter and antimatter. 
As we mentioned at the beginning of this section, the antimatter domain should be very 
far at ^ > Gpc. Moreover, there is another problem with this mechanism, namely domain 
walls between matter and antimatter domains could destroy the observed homogeneity and 
isotropy of the universe [60]. To avoid the problem a mechanism of the wall destruction 
is necessary. 

3 Stochastic or dynamical. If a complex scalar field \ was displaced from its equilibrium 
point in the potential, e.g. by quantum fluctuations at inflation, and did not relaxed down 
to equilibrium at baryogenesis, it would create CP-violation proportional to the amplitude 
of the field but without problems of spontaneous CP-violation. Later, after baryogenesis 
was over, % would relax down to zero. So domain walls do not appear. Inhomogeneous 
P(x) with domains of matter and antimatter can be created with such CP-violation. Their 
size depends upon the details of the scenario. 

There is a long but probably incomplete list of different scenarios of baryogenesis (BG): 

1. Heavy particle decays [54] . 

2. Electroweak BG [HI]. Too weak in MSM but may work with TeV gravity. 

3. Baryo-through-leptogenesis [tj2"] . 

4. SUSY condensate BG [63]. 

5. Spontaneous BG [64] . 
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6. BG by PBH evaporation [65] . 

7. Space separation of B and B at astronomically large distances jHS], which is probably 
not effective. However anti-baryons might be removed from our into higher dimensions [HZ] 
or predominantly accumulated inside quark nuggets [68J . 

7. BG due to CPT violation [71]. 

In all these scenarios new physics beyond minimal standard model is necessary In what 
follows we will very briefly describe some of these scenarios. More details can be found in 
the reviews |72j . 

Heavy particle decay BG is naturally realized in grand unification theories, GUTs, 
where gauge bosons X with mass around 10 16 — 10 15 GeV are present. These bosons can 
decay e.g. into, qq and ql pairs where baryon number is evidently not conserved. Due to 
large mass of X the deviation from equilibrium could be significant, CP- violation might be 
sufficiently large (we know nothing about it) and the mechanism could be efficient enough 
to generate the observed asymmetry. The problem with GUTs is that the temperatures 
of the GUT scale might not be reachable after inflation. On the other hand, baryogenesis 
might proceed with under-abundant X-bosons created out of equilibrium. 

Particles and antiparticles can have different decay rates into charge conjugated chan- 
nels if C and CP are broken, while the total widths are equal due to CPT invariance. If 
only C is broken, but CP is not, then partial widths, summed over spins, are the same 
because CP-invariance implies: 

T(X^f,a) = (X^f,-a). (199) 

If both C and CP are broken, partial widths may be different, but the effect takes place 
happen in higher orders of perturbation theory. In lowest order the amplitudes of charged 
conjugated processes must be equal, A = A*, because of hermicity of Lagrangian. The 
same would be also true for higher order contributions if they were real. An imaginary 
part is generated by re-scattering in the final state (with non-conservation of B or L), as 
can be seen from the S-matrix unitarity condition: 

i(T if - 4) = - T inTl } = -£ T} n T n} (200) 

n n 

Let us consider an example of X-boson decays into the channels: 

X -> qq, X -»• ql, 

X^qq, X^ql. (201) 
and assume that the partial widths are different due to C and CP violation: 

r x ^ qq = (i + A q )v q , r M - = (i - A,)r f , 

T^ qq = (1 - A q )T q , Fx^i = (1 + A l )F l . (202) 

Here r ~ a and A ~ a, where a ~ 1/50 is the fine structure constant at GUT scale. 
The asymmetry is proportional to /3 ~ (2/3) (2 A 9 — A;). Its magnitude can be roughly 
estimated as 

n Sf AT m , , 

P ~ -r - p 203 

/ r mpi 
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Small numerical coefficients omitted here would diminish the result. For example, the 
subsequent entropy dilution by about 1/100 is not included. For successful lepto/baryo- 
genesis the mass of the decaying particle should be larger than 10 10 GeV, or m P i <C 10 19 
GeV). 

Problem 21. How the charge asymmetry generated in heavy particle decay vanishes 
in equilibrium? It is stated in the literature that the inverse decay does the job. However 
one can see that it is not so because using CPT, one finds: 

v m ^ x = (i + A,)r 9 , = (1 - A,)r,, 

v qq ^ x = (i - A,)r 9 , r ql _ x = (l + A,)r,. (204) 

Thus direct and inverse decays produce the same sign of baryon asymmetry! 

Electroweak baryogenesis is very attractive because all the necessary ingredients are 
present in the minimal standard model. CP is known to be broken, but, as we have 
seen above, very weakly. Baryonic number is non-conserved because of nonabelian chiral 
anomaly. At zero T baryon nonconservation is exponentially suppressed as exp(— 2ir / a) [73] , 
because of barrier penetration between different vacua. However, it is argued that at 
high T it is possible to go over the barrier, by formation of classical field configuration, 
sphalerons. We do not know how to calculate the probability of production of large coher- 
ent field configurations in elementary particle collision but lattice simulations show that 
their production might be efficient and sphalerons could have thermal equilibrium abun- 
dance. Unfortunately the deviation from equilibrium of massive particles at EW scale 
is tiny and the first order electroweak phase transitions seems to be excluded because of 
heavy Higgs boson. So electroweak baryogenesis, though very attractive is not efficient 
enough, but it may be operative with TeV gravity. 

Baryo-through-leptogenesis is probably the most popular mechanism today. It started 
from creation of lepton asymmetry by L-nonconserving decays of heavy, m ~ 10 10 GeV, 
Majorana neutrino, analogously to GUT, and subsequent transformation of the lepton 
asymmetry into baryonic asymmetry by CP symmetric and B non- conserving and (B — L) 
conserving electroweak processes. The mass matrix of three flavour light and heavy Ma- 
jorana neutrinos has 6 independent phases, three in the sector of light neutrinos and 3 in 
heavy ones. They are unknown and allowed to be of order unity. 

Primordial black hole evaporation does not demand B-nonconservation at particle 
physics level for generation of the baryon asymmetry. Of course, thermal evaporation 
cannot create any charge asymmetry. However the spectrum is not exactly black but 
is modified due to propagation of the produced particles in gravitational field of BH. 
Moreover, an interaction among the produced particles is essential. Let us assume that a 
meson A is created at the horizon and decays as: 

A^ H + L and A -> H + L, (205) 

where H and L are a heavy and light baryons, e.g. t and u quarks, respectively. Due to 
CP- violation the branching ratios of these decays may be different. Back-capture of H 
by gravitational field of the black hole is larger than that of L. Thus some net baryon 
asymmetry in the external world could be created. If the cosmological energy density 
of black holes at production was small, Pbh/ ptot = e < 1, then at red-shift z — 1/e, 
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with respect to the production moment, the non-relativistic BHs would dominate. Their 
evaporation could provide the necessary baryon asymmetry and reheat the universe. 

For a numerical estimate let us present some simple formulae, omitting numerical 
factors of order 1 — 10. The black hole temperature is essentially given by the only 
available parameter with dimension of length, i.e. by its gravitational radius: 

T BH ~ l/r 9 ~ m 2 Pl /M BH . (206) 
The luminosity of the body with temperature Tbh and radius r g is: 

Lbh ~ T A r 2 g ~ m%jM 2 BH . (207) 

Correspondingly the BH life-time is equal to: 

tbh ~ M% H /m% v (208) 

For example, if Mbh = 10 15 g, its life-time is equal to the universe age, tbh ~ tu ~ 10 10 
years. For our case much lighter BHs are needed. Let us assume that primordial BHs 
were formed at Tbh = 10 14 GeV and their mass was equal to the mass inside the cosmo- 
logical horizon at that moment, Mbh — m 'pit ~ 10 4 g- The life-time of such BHs would 
be Tbh ~ 10~ 16 sec, which corresponds to cosmological temperature T ~ 10 5 GeV and 
red-shift from the moment when horizon mass was equal to Mbh, was about 10 9 . In 
other words, if the mass fraction of BHs at production was 10 -9 , then at the moment of 
their evaporation they would dominate the cosmological energy density and could create 
observed baryon asymmetry even if the fraction of the baryon number density was small 
in comparison with the total number density of the evaporated particles. 

Spontaneous baryogenesis may operate in thermal equilibrium. Explicit CP-violation 
is not obligatory. It is assumed that a global [/(l)-symmetry associated with baryonic 
number is spontaneously broken. The Higgs-like scalar boson acquires non-zero vacuum 
expectation value and its phase becomes massless Goldstone boson, = n exp(i0). In the 
broken phase the Lagrangian can be written as: 

C = v 2 (d9) 2 + dflil - V{9) + iQ^QiL^L + (gnQL + h.c). (209) 

In the case of homogeneous 9(t) the second term looks like chemical potential, On^. How- 
ever, in reality it is not true, because chemical potential is introduced into Hamiltonian 
but for derivative coupling £ ^ Ti. 

If the potential V{6) = 0, i.e. in purely Goldstone case, we can integrate the equation 
of motion: 

2n 2 d 2 6 = -dtf (210) 

and obtain: 

An B = -r] 2 A9, (211) 

i.e. non-zero baryon asymmetry in thermal equilibrium and without explicit CP- violation. 
The latter is created by initial 9 ^ 0. 
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In realistic situation 9 is small (because inflation kills all motion) and the pseudogold- 
stone case, i.e. non-zero V(9) could be more efficient. Now the equation of motion for 9 
takes the form 

r ] 2 9 + 3H9 + V'(9)=d fl jZ, (212) 

where V(9) ~ m 2 if [—1 + (9 — it) 2 ] and = ^7 M ^. Initially 9 is uniform in [0, 2tt] and 
after inflation it started to oscillate around minimum. 

The second necessary equation is that for the quantum baryonic Dirac field: 

(id + m) rj) = -grjl + {d„9)j^ (213) 

The solution to this equation can be found in one- loop approximation for ip(0) in external 
classical field 9. Then this solution, ipty = F(9), should be substituted into eq. (12121) . 
In this way a closed equation for 9(t) can be obtained. The solution oscillates with 
alternating baryonic number giving the net result for the baryon number density 

n B ~ v 2 T AB (A9) 3 . (214) 

The SUSY baryonic condensate scenario will be discussed in more detail here because 
with simple modification it allows for creation of astronomically significant antimatter [69] . 

The basic features of this scenario are the following. SUSY predicts existence of scalars 
with non-zero baryonic number. Such bosons may condense along flat directions of the 
potential: 

tA(x) = A|x| 4 (l-cos4£), (215) 

where \ = \x\ exp(z'0). In SUSY models with high energy scale the baryonic number is 
naturally non-conserved. It is reflected by the non-sphericity of potential ( 12 15ft . Due to 
infrared instability of massless (m <C H) fields in de Sitter space-time, x can travel away 
from zero along the flat directions, 9 = 0, tt/2, it, 3tt/2. We can also add a mass term to 
the potential: 

U m ( X ) = m 2 \ X \ 2 [1 - cos(2# + 2a)] , (216) 

where m = \m\e a . If a ^ 0, then C and CP are explicitly broken, though it is not necessary 
for baryogenesis, 

"Initially" (as a result of inflation) x was pushed away from origin and when inflation 
was over it started to evolve down to equilibrium point, x = 0? according to the equation 
of the Newtonian mechanics: 

X + 3H X + U'(x) = 0- (217) 

The baryonic number of X'- 

B x = 9\x\ 2 (218) 

is analogous to mechanical angular momentum. Using this mechanical analogy, and having 
the picture of potential U (x) is easy to visualise the solution of the equation of motion 
without explicitly solving it. 
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The baryonic number of x is accumulated in its "rotational" motion, induced by 
quantum fluctuations in orthogonal to valley direction. When \ decays its baryonic 
charge is transferred to that of quarks through B-conserving processes. If the mass term 
is absent or symmetric with respect to the phase rotation of x this scenario leads to 
globally charge symmetric universe. The domain size Ib is determined by the size of the 
region with a definite sign of 9. Usually Ib would be too small if no special efforts are 



If m 7^ 0, the angular momentum, B, is generated by a different direction of the mass 
valley at low x- If CP-odd phase a is small but non-vanishing, both baryonic and an- 
tibaryonic regions are possible with dominance of one of them. In this case matter and 
antimatter domain may exist but globally B ^ 0. 



Now let us modify the model by adding general renormalizable coupling of x to inflaton 
field $: 



where $i is the value of $ which it passed during inflation, not too long before its end. 
It is a free adjustable parameter. 

Because of this coupling the gates to the valley would be open only for a short time 
when $ was close to $i. So the probability for x to reach a large value would be small. 
As a result we will have the following picture of the universe. The bulk of space would 
have normal homogeneous baryon asymmetry, (3 = 6 • 10 -10 , with small bubbles having 
large f3 ~ 1. In the simplest version of the scenario the high B regions should be almost 
symmetric with respect to baryons and antibaryons. 

The mass spectrum of such baryon rich bubbles is practically model independent (it 
is determined by inflation) and has simple log-normal form: 



Such object could make primordial black holes, quasars, disperse clouds of antimatter, 
and unusual stars and anti-stars, all not too far from us in the Galaxy. Phenomenolog- 
ical implications of this mechanism of antimatter creation and observational bounds are 
discussed in ref. [70]. If such mechanism is realized in nature the attempts for search of 
cosmic antimatter have non-zero chances to be successful. 

It is worth noting that primordial nucleosynthesis in high B domains proceeded with 
large ratio ris/n^. So the outcome of light and heavier element abundances could be 
much different from the predictions of the standard BBN. Thus the regions in the sky 
with abnormal chemistry would be first candidates to search for cosmic antimatter through 
100 MeV photons or through the positron annihilation line. 



done. 



A $ | X | 2 (*-*i) 2 , 



(219) 
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